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The implications of restricted conformal invariance under conformal transformations 
preserving a plane boundary are discussed for general dimensions d. Calculations of the 
universal function of a conformal invariant £ which appears in the two point function of 
scalar operators in conformally invariant theories with a plane boundary are undertaken to 
first order in the e = 4 — d expansion for the the operator </> 2 in </> 4 theory. The form for the 
associated functions of £ for the two point functions for the basic field (j) a and the auxiliary 
field A in the the iV — > oo limit of the O(N) non linear sigma model for any d in the range 
2 < d < 4 are also rederived. These results are obtained by integrating the two point 
functions over planes parallel to the boundary, defining a restricted two point function 
which may be obtained more simply. Assuming conformal invariance this transformation 
can be inverted to recover the full two point function. Consistency of the results is checked 
by considering the limit d — > 4 and also by analysis of the operator product expansions 
for cf) a (f)P and AA. Using this method the form of the two point function for the energy 
momentum tensor in the conformal O(N) model with a plane boundary is also found. 
General results for the sum of the contributions of all derivative operators appearing in 
the operator product expansion, and also in a corresponding boundary operator expansion, 
to the two point functions are also derived making essential use of conformal invariance. 



1 Introduction 



In more than two dimensions it is not generally possible to construct explicitly non 
trivial conformal field theories, with a detailed knowledge of the spectrum of spins and 
scale dimensions of all operators in the theory and further the coefficients appearing in 
operator product expansions for each pair of operators, at least to the same degree as in 
two dimensions. Nevertheless for a very large class of quantum field theories scale invari- 
ance at possible critical points may also be extended to invariance under the full conformal 
group [1] which implies significant restrictions on the form of multi-point correlation func- 
tions [2]. In particular the functional form of two and three point functions is effectively 
unique assuming conformal covariance since to construct conformal invariants, and hence 
for arbitrary functions to be present, four or more points are required (for operators with 
spin there may be more than one linearly independent conformally covariant form for the 
three point function [3]). From an experimental viewpoint in a statistical physics context 
only two point functions are mostly relevant and in this case conformal invariance gives 
little more than just scale invariance. On the other hand for calculating critical exponents 
in the 1/N expansion using conformal invariance has proved essential in obtaining results 
to 0(1/N 2 ) and 0(l/iV 3 ) [4]. 

For statistical mechanical problems involving a boundary then at a critical point there 
are new critical exponents, expressing the behaviour of physical quantities near or on the 
boundary, which are unrelated (at least in any simple general fashion) to bulk critical 
exponents [5]. Also for any particular bulk critical point there are a variety of possible 
boundary conditions with differing surface exponents. As Cardy first showed [6,1] there is 
still a residual conformal group consisting of conformal transformations leaving the bound- 
ary invariant. For a plane boundary in d Euclidean dimensions the restricted conformal 
group is then 0(d — 1,1). In this case the two point function, involving operators at x, x', 
depends on functions of a single conformal invariant £(x,x'). These functions depend on 
the particular theory and associated boundary conditions, or rather on their correspond- 
ing universality class. However conformal invariance is a significantly stronger requirement 
than scale invariance as far as potential experimental implications are concerned [7]. If 
we define coordinates = (y, x), with y measuring the perpendicular distance from the 
boundary, then scale invariance by itself, with conventional translational and rotational 
symmetries, only restricts the two point function to depend on functions of the two scale 
invariant variables s 2 /y 2 , s 2 /y' 2 where s = x — x'. 

For scalar fields there is a single function of £ in the associated two point function 
but for fields with spin there may be several. As a particular illustration we consider the 
energy momentum tensor T MZ ,, which is traceless in the conformal limit, and the two point 
function then contains three possible invariant functions. However the conservation equa- 
tion dfj,Tfj, v = provides two first order linear differential equations linking these functions. 
In two dimensions the number of invariant functions is reduced to two and in this case the 
differential equations have a unique solution, satisfying appropriate boundary conditions, 
in accord with known results [8]. Under the restricted conformal group it is also possible 
to form non zero two point functions for fields of differing spin and scale dimensions, such 
as for and a scalar field although in this case the functional dependence is entirely 
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determined. 

In this paper we investigate conditions on the form of two point functions for operators 
in conformally invariant theories with a plane boundary. Previously we discussed such 
two point functions involving the energy momentum tensor and derived the necessary 
conditions arising from the conservation equations for T^ v by considering a perpendicular 
configuration where the two points were restricted to be on a perpendicular to the boundary 
[9] (referred to subsequently as I). Here we construct a conformally covariant form for the 
general two point function {T^ v {x)T pcr {x')) for arbitrary x,x' by constructing in terms of 
the invariant scalar £(x,x') vectors X p oc <9^£, X' a oc d' a ^ at x,x' respectively. 

In general a two point function in the presence of a boundary may be written, 
expressing invariance under translations parallel to the boundary, as (0(x)0'(x')) = 
G(y, j/',x — x'). In our discussions it is useful to consider a transformation of G obtained 
by integrating over planes parallel to the boundary, f d d_1 x G(y, y', x) = G\\ (y, y'), which 
defines what is sometimes referred to as the parallel susceptibility. For a conformally in- 
variant theory, when (0(x)0'(x')) is expressible in terms of a function #(£), this procedure, 
which we refer to as parallel integration, defines a new function g(p), p = (y — y') /Ayy' . 
Crucially we are able to show that the transformation g ^ g is invertible with g expressible 
in terms of an integral over g, for d = 3 the result is simply #(£) = — g'(£)/ir. This then 
allows for significant simplifications in calculations since finding G\\ (y, y') is sufficient to 
determine the full two point function (O(x)O'(x')). 

These methods are illustrated by application to the O(N) cr-model in the large N 
limit. This is in the same universality class as an N component scalar field theory with 
a renormalisable 4> 4 interaction at the non Gaussian Wilson fixed point. We are able to 
recover some results obtained some time ago for the arbitrary function of £ associated with 
the two point functions of the basic N component fields <j> and also the auxiliary field A 
which is a scalar under O(N) [10,11]. 

We also consider constraints on the invariant function #(£) arising from the operator 
product expansion in which the operators with lowest dimension and non vanishing one 
point functions are relevant for the limit £ — > 0. With conformal invariance only scalar 
operators, or their derivatives, in the operator product expansion of O and O' can con- 
tribute to their two point function in the presence of a boundary. We derive an explicit 
form for #(£), in terms of hypergeometric functions, resulting from all derivative operators 
formed from any particular scalar operator occurring in the operator product expansion. 
The derivative terms in the operator product expansion are determined by the requirement 
that they should reproduce exactly the appropriate full three point function for the con- 
formal field theory without boundary. These results are applied to two point amplitudes 
involving T^ v and the expressions obtained automatically satisfy the required conservation 
equations. 

A similar boundary operator expansion, where a bulk operator 0(x) is expanded in 
terms of boundary operators 0{~x) [12], is also investigated. This constrains the function 
#(£) in the limit £ — > oo. Again, using conformal invariance, we are able to sum up 
explicitly the results for all derivative operators formed from a given boundary operator 
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which in general also involves hypergeometric functions. In this case the essential input 
determining the derivative terms for a particular boundary operator O appearing in the 
expansion of O is the form of the two point function for O and also O at arbitrary points. 
The boundary operator expansion of the energy momentum tensor defines a boundary 
scalar operator T which is given by the non singular limit of Tj_±(x) as x — > (0, x) and 
therefore has dimension d [8] . However for d > 2 it is necessary to also consider boundary 
operators %j, symmetric traceless tensors whose dimension is not constrained by general 
principles (this again reflects the non uniqueness of (T flv ,(x)T ap (x')) for d > 2). 

In detail in section 2 we consider the general conditions stemming from conformal 
invariance for two point functions with a boundary. In particular we describe the notion 
of quasi-primary operators which have simple properties under conformal transformations 
and analyse in detail the forms of the two point amplitudes for vector operators V M and also 
T pu in the presence of a boundary. Using the vectors X^, X' a it is straightforward to write 
down conformally covariant expressions. In section 3 we consider the simplest conformally 
invariant scalar field theories defined by free fields and also with a quartic interaction at the 
non Gaussian fixed point present in the e = 4 — d expansion. In particular we calculate the 
two point function for the operator </> 2 to first order in e. In section 4 we analyse the O(N) 
model in the large N limit and describe in detail the transformation g <-> g mentioned 
above. Section 5 contains more calculations for the two point functions of U M or T^ v in 
the O(N) model to leading order in 1/N. The result for the two point function of the 
energy momentum tensor is non trivial, involving generalised hypergeometric functions. 
Sections 6 and 7 respectively contain the details of our discussions of the consequences of 
the operator product and boundary operator expansions with applications to the results 
obtained in the O(N) model. Various calculational details are relegated to five appendices. 

2 Conformal Invariance with Plane Boundaries 

In flat d-dimensional Euclidean space with coordinates conformal transfor- 

mation g is defined by preservation of the line element up to a local scale factor 

x„ -> x\(x) , dx^dx^ = 3 (x)" 2 dx M dx M d d x 9 = Q 9 (x)- d d d x . (2.1) 

Such transformations define a group, (x 91 ) 92 = x 9291 , which is isomorphic to 0(d + 1, 1). 
For any conformal transformation we may define a local orthogonal matrix belonging to 
O(d) by 

n 9 ^(x) = n 9 (x)^ , (2.2) 

satisfying 1Z 92 (x 9l )TZ 91 (x) = TZ 929l (x). For d > 2 arbitrary conformal transformations can 
be generated by combining translations and rotations, for which Q 9 = 1, with inversions 
through the origin, i, taking x^ — > x^/x 2 , for which 

^VW = VW^f-2^, n l (x)=x 2 . (2.3) 

Under conformal transformations I liv {x — x') — > TZ fl0 t(x)TZ 1/ /3(x')I a p(x — x') so that it 
transforms as a vector at x and also at x'. In consequence I liv {x — x') may be regarded 



3 



as representing a form of parallel transport for conformal transformations and it plays a 
crucial role in the construction of two point functions for vector and other tensor fields. 

In a conformal field theory fields O l (x) which form the vector space on which some 
finite dimensional irreducible representation of 0(d) acts are defined as quasi-primary fields 
if they transform under the full conformal group according to [13] 

0\x) -> O ai (x 9 ) = ^(x^D^TZ^x^O^x) , (2.4) 

where rj is the scale dimension of the fields and, for any R e 0(d), D(R) is the corre- 
sponding matrix in this representation. Since the transformation rule in (2.4) depends on 
x, derivatives of quasi-primary fields are not in general quasi-primary but transform with 
extra inhomogeneous terms. However for V M a vector field of dimension d — 1 or T^ v a 
symmetric traceless tensor field of dimension d it is straightforward to show [3] that d^V^ 
and d^T^ are quasi-primary fields. Hence it follows that in a conformal field theory a con- 
served vector current and the energy momentum tensor necessarily have dimensions d — 1 
and d respectively. For the full conformal group there are no invariants which can be con- 
structed from two or three points so that the two or three point functions of quasi-primary 
fields are essentially uniquely determined, with no arbitrary functions present. 

For a flat Euclidean space with a plane boundary, as was shown by Cardy [6] , a non 
trivial subgroup of conformal transformations still remains. If we let x M = (y, x) and define 
a plane boundary by y = 0, so that y is the perpendicular distance from x to the boundary, 
then it is necessary to restrict the conformal group to those transformations leaving y = 
invariant. This subgroup is then generated by d— 1 dimensional translations and 0(d — 1) 
rotations acting on x together with the inversion x^ — > x^/x 2 again and forms the group 
0(d, 1). Under such transformations for two points a; M , x'^ it is easy to see that 

/ _ /\2 _ (x — x ) 2 y , y' , . 

[x x ' ~* n(x)n(x') ' y ~* n(x) ' y ~* Q( X ') ■ 1 ' 

Hence we may construct invariants 

= (X-X'f 2 = (X-X'f = 

? 4yy' ' (x - x') 2 + 4yy> £ + 1 " 1 ' ' 

where < £ < oo and < v < 1. 

For a one point function of an operator O in the neighbourhood of a plane boundary 
conformal invariance under transformations as in (2.4) implies that this can only be non 
zero for quasi-primary scalar fields, belonging to the singlet representation of 0(d), when 
we can write 

w = • < 2 - 7 > 

For the two point function of quasi-primary operators the existence of the conformal 
invariants in (2.5) implies that there may be an arbitrary function present. In addition the 
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two point function may be non zero for operators of differing spins and scale dimension, 
unlike the case for conformal invariance without a boundary. For two scalar fields we may 
write in general 

= { 2y)H2y')^ /l2(0 = \x-x')^ F ^ ' ^MO = FMv). (2.8) 

The functions f 12(C) or F 12 (v) are constrained by the operator product expansion (OPE). 
If for s = x — x' — > 

Oi(x)0 2 (x') Mx') , (0 3 (x')) = -4^— , (2.9) 

then for £ ~ — > there is a contribution to the functions f 12(C) or Fi 2 (v) in (2.8) of the 
form 

/12(C) ~ C^AsC^ 112 -^ , F 12 (v) ~ ds 3 ^ 1 "^ 3 . (2.10) 

Later, in section 6, we determine the entire contribution of non leading derivative terms 
in the OPE. 

There are additional constraints on the functions f 12(C) or Fi 2 (v) arising from a 
boundary operator expansion (BOE). For conformal field theories with boundary it is nat- 
ural to define a basis of boundary operators 0„(x), which transform irreducibly under 
0(d — 1) and have a well defined scale dimension rj n . The precise set of such operators 
depends of course on the particular conformal theory and also the precise boundary con- 
ditions, compatible with conformal invariance imposed. The BOE has the form 

where in (2.7) Aq = Bo 1 . The two point function of a bulk operator and a boundary 
operator is determined up to an overall constant 

(0(x)6„(x-)) = (2y) f_°f ( -. 2) ,„ ■ V = fe,x-x'). (2.12) 

Combining (2.11) and (2.12) it is easy to find that in the limit C — ► 00 or v — > 1 

/ 12 (0 ~ B 0i6 B Sc- f ' , F 12 (v) ~ B 0i6 B 02 6 (\ - v 2 f~^ , (2.13) 

where fj is the scale dimension of the leading boundary operator O which gives a non 
zero contribution when using (2.11) and (2.12) to take the limit y' — > 0. Clearly there are 
consistency conditions arising from considering the alternative limit y — > 0. For a statistical 
mechanical system at a critical point the invariant functions /(£) or F(v) appearing in two 



5 



point functions are universal, depending only on the universality class and the boundary 
conditions.* 

For non scalar quasi-primary fields the expressions for the two point function in the 
neighbourhood of a plane boundary are more complicated to analyse. In I a conformal 
transformation was used to restrict x M = (y, 0) and x'^ = (y',0), on a perpendicular 
to the boundary, where the remaining invariance or little group is reduced to 0(d — 1) 
rotations. However we now describe an alternative approach which gives equivalent results. 
This is based on defining vectors X^, X'^, with zero scale dimension, under restricted 
conformal transformations preserving the boundary at x, x' respectively, so that X^ — > 
TZ fia (x)X ce , X — > 7Z^ a (x')X'i3. Since £ in (2.6) is a scalar these are explicitly given by 

= yjd^ = ~ n i*) » x 'i* = y'j 9 '^ = v {--^- s » - n i*) > ^ = (1,0), 

(2.14) 

which satisfy 

X^X^ = X ^X M = 1 , X^ = I^ v (s)X v . (2-15) 
As y -> X^ -> -ra^ and for y' -> A^ -> -I^ v (s)n u . 

We may now easily construct invariant forms for two point amplitudes for tensor fields. 
For a vector field of dimension d — 1 then 

(^(aOKCaO) = j^p^W^ + Wr)) . (2.16) 



To impose current conservation, d^V^ = 0, we use 



2y' 



2W 



^(^pW*))= . ^((^t x m) =-(d-i) (g2) 

<9 M A'„ = -^^(-/^(a) + , dfj,F(v) = -^X^F'iv) , 

s s 



(2.17) 



to obtain 

v^-{C + D) = {d-l)D. (2.18) 

Away from the boundary, for v — > 0, .D — > while C is a constant determined by the bulk 
conformal theory. Similarly for T^,, the traceless energy momentum tensor of dimension 
d and O a scalar field of dimension r\ we may find expressions for the mixed two point 
functions 

( c 2i/) d ~ 1 ~ r i 
(V,(x)0(x')) = { J) X,C vo {v), 

[ \ (2-19) 
(2v') d ~ 71 / 1 \ 

(T^(x)0(x')} = K -f^[X,X v - -S^C T o(v) . 



* Note that the results (2.8) give for the limits of the two point function perpendicular and parallel 
to the boundary (Oi(y, 0)O 2 (y' , 0)) oc l/y v± as y -> oo and (Oi(y, x)0 2 (y, x')) oc l/|spl as |s| -> oo 
where the surface critical exponents ri± = rji + r), r]\\ = 2f) [10]. 
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Current conservation and also d^T^ = 0, using (2.17) and 

a,(^(lA-^))=-(rf-l)^T^, (2-20) 

requires 

v^-Cvo = (d- l)C vo , v^-Cto = dC T o , (2.21) 
dv dv 

which gives a unique functional form 

C vo (v) = c V ov d - 1 , Cto(v) = c TO v d . (2.22) 

The coefficients cyo and cto are determined by Ward identities. The energy momentum 
tensor T Mi , may be defined by the response to arbitrary infinitesimal reparameterisations 
Xfj, — > x^ + a fl (x) [6,14] and from conformal invariance it is possible to derive the relations 
(see eqs. (E.7,8) in ref. [3]) 

(T^(x)0(x')) = dC(0(x'))8 d (s), 
d»{T^{x)0{x')) = (1 + C) (0(x'))d„5 d (s) - (d„0{x'))5 d {s) , ^ 

where C is undetermined (reflecting the arbitrariness of {T^ u {x)0{x')) when regularised 
up to terms oc 5^5 (x — x')(2y')~ ri ). The expression obtained in (2.19,22) can be shown 
to be compatible with (2.23) by considering the expansion for s — > which has the form 



' d 1 I f 1 



- -^(n^s„ + n u s^- 8 flu n-s+ (d-2)n-s ^ 

With a suitable regularisation of the terms 0(s~ d ) we may find * 
v d /„ „ 1 , \ d - 1 „ 1 „ ^, „ 1 



(2.24) 
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^ - \*») = - (^F a ^ d(x - x>) + c wr dJd{x ~ x>) ' (2 - 25) 

for c arbitrary and Sa = 2n^ d /T(|<f). Assuming the result (2.25) the regularised trace of 
(2.24) is then also dc5 d (x — x')/{2y) d . With these formulae it is easy to check that the 
identities (2.23) are satisfied, using (2.7), if CA G = ctqc and [8,7] 



* In the spirit of differential regularisation [15] we may obtain a well defined distribution for the 

0(s~ d ) terms -\ ( 1t2±L _ I § \ = (dnd v — -S^vd 2 } } - + c5u, v 5 d (s), where c is an arbitrary 

s a \ s z a J d(d — 2) V d J s a z 

constant. The divergence and trace of this expression may be found using — d 2 s~ d+2 = (d — 2)SdS d (s). 

The trace of the 0(s~ d+1 ) terms in (2.24) is zero while the divergence is unambiguous being oc S d (s). 
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The two point function of the energy momentum tensor itself can also be written in 
a similar fashion in terms of a basis of conformally covariant symmetric traceless tensors 
at x,x' in the general form 

(T^(x)T ap (x')) = -^y^u,a P (s)C(v) + (X„X V - ^) (x' a X' p - \^ p )Mv) 

+ [x^X'J up {s) + /i <-> v,a <-> p (2.27) 

4 / / 4 4 \ "I 

-d^X'tjX'p - -dapX^Xv + -pS^Sap) B(v) > , 

where T^^p represents inversion on symmetric traceless tensors 

1^a P {s) = \ (l^{s)I vp {s) + I^ p {s)h a {s)) - - 8^8 ap . (2.28) 

If the two point function is defined on all M d , with no boundary, then A, B are absent while 
C is a constant Ct and the conservation equation for T^ v , giving dp i {T lxv {x)T ap (x')) = 0, 
is satisfied by virtue of 

^((^W^O- (2-29) 
In the general case the conservation equation, using (2.29) with (2.17), (2.20) as well as 

2y> v (2 ' 30) 
~{{2-d 2 )I V(T {s) + {d-2)X v X' (T ), 



(s 2 ) d+1 d 
leads to the conditions 

(v4- ~ d) (C + 2B) = - -(A + AB) - dC , (2.31a) 
\ av / a 

(v— - d) ((d - 1)A + 2(d - 2)B) =2A- 2(d 2 - A)B . (2.316) 

These equations still leave the A, B, C undetermined up to an arbitrary function of v. 
Away from the boundary, or as — > 0, the terms involving X^X' a should disappear, 
hence in this limit A, B — > 0, while C — > Ct which is the constant determining the scale 
of energy momentum tensor two point function which has a simple form in the associated 
conformal field theory without boundary, 

Ct 

\I'p,v{p^)T cr p( K X ))no boundary 1 ^p.irpf^) • (2.32) 

\ S ) 

On the boundary if we assume Tu(0, x) = 0, with % denoting components tangential to the 
boundary, then we must require C(l) +2.6(1) = 0. Since (T MI/ ) = with a plane boundary 
there are no Ward identity relations for this two point function. 
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When d = 2 there are only two linearly independent tensors instead of the three in 
(2.27) so that only A, C + 2B are relevant. The differential equations then have a unique 
solution obeying the boundary conditions 

A(v)=4C T v 4 , C(v) +2B(v) =C T (l-v 4 ) . (2.33) 

In two dimensions with complex coordinates z = x + iy, z = x — iy, 5 ZZ = \ then 

X z X' z = \l zz {s) = -- Z -^, Xz X'- z = - Z -^, l z - z (s) = 0, (2.34) 
and since v 2 = \z — z'\ 2 /\z — z'\ 2 we thereby obtain the universal form [1] 

(T zz (x)T zz (z')) = , {T zz {z)T,- z {z>)) = . (2.35) 

Up to a suitable normalisation factor Ct is of course the Virasoro central charge. 

For general d the expression (2.27) for the energy momentum tensor two point function 
may be related to the results in I, where free field expressions were given for arbitrary d 
and also the functional dependence on v calculated to first order in the e expansion from 
d = 4 for the non Gaussian fixed point in (p 4 field theory, by restricting (2.27) to the 
perpendicular configuration, when for y > y' ', = —X'^ = n^, v = (y — y') / (y + y'), the 
general form can be written as 

(Tij(v,0)Tu(l/,0)) = 7- 4j 7(^, 

{y " ' (2.36) 
(Tijiy, 0)T ke (y',0)) = _ (S(v) 5 l3 5 u + e(v)(6 ik 6 je + 6 a S jk )) ■ 

All other non zero components may be obtained trivially from (2.36) using T^ v = T V[X and 
= 0, for instance 

(T n (y, OjTuCy', 0)) = _^ a(v) , a = (d - l)((d - 1)5 + 2e) . (2.37) 



It is straightforward to show that 

a= d -^ L {(d-l)(A + AB)+dC) , 1 =-B-\C, e = ±C, 5 = -^(A + AB - dC) . 



(2.38) 

Eqs. (2.31a,b) then translates into the results given in I, for instance (2.31b) becomes 
more simply 

(v-^-d^ja(v) = 2(d-l)^(v). (2.39) 



On the boundary itself it is easy to see that 

1 

(s~*) 



(Tn^xjTn^x')) = 7^a(l) . (2.40) 



3 Scalar Field Theory 



The simplest conformally invariant field theory is that given by free massless scalar 
fields 4> a {x), a = 1 ... N (for later convenience we assume N components). To ensure com- 
patibility with conformal invariance we may impose either Neumann or Dirichlet boundary 
conditions which for a plane boundary at y = requires <9i(/> a (0,x) = or (/> a (0,x) = 
respectively. If 

(0 a (z)0V)> = S^G^x, x') , G^x, x') = ^F^v) , (3.1) 
then in the free field case it is easy to see, using the method of images, that 

7to = §d-l, F 4> {v)=A(l±v d - 2 ), (3.2) 
where the upper/lower signs correspond to Neumann/Dirichlet boundary conditions and 

A= w J w d - Sd= Y(id)- (3 - 3) 

For the composite operator 4> 2 with dimension 77^2 then in general 

^\x)<?{x>)) = (<p 2 (x))(<p 2 (x')} + G^(x,x'), 

G^(x,x') = -1-F^iv), ((P 2 (x)) = (3 ' 4) 

and in the free field case 

V(j) 2 =2 V4> = d-2, F^(v) = 2NF (f) (v) 2 , A^=±NA. (3.5) 



A more interesting case is the conformal field theory realised at the non Gaussian 
fixed point in the theory with interaction ^g((j) 2 ) 2 for which critical exponents and other 
universal quantities can be calculated in the e = 4 — d expansion. With minimal subtraction 
at the fixed point g*/lQn 2 = 3e/(N + 8) + 0(e 2 ). For the (cfxfi) two point function the 0(e) 
corrections to the free field results for the universal function F^v) have been calculated 
by Gompper and Wagner [16] and also in I*. Since g* = 0(e) it is sufficient to first order 
to evaluate the Feynman integrals just for d = 4. The results obtained at one loop give 
774, = \d - 1 + 0(e 2 ) and 

F ^ (1, = Tl^ s As(vH n 1 -^±Ml-v<)) . (3.6) 



* In I the results in (C.1,2), and also (C.3), for G(y,y') should have an additional ± sign. Other 
misprints in I are that the factor multiplying F" in (C.5) should be v 2 (l — v 2 ) 2 , in (C.8a) the tensor 
structure should be 5ik5j£ + SaSjk and in the last line of (C.8b) the second factor should be x 2 + z 2 — y' 2 , 
in (2.6) in the transformation equation for T^ v replace Cl(x) by £l(x) d , in (2.26) the correct equation is 
Aij = r)5ij/(d— l)Sd and in (2.29) in the formula for a(l) 2d — > 2 d . Further in (A. 2), and also subsequently 
on the same page, e^ie^j — > e fJ, ie u j. 
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In the Dirichlet case it is easy to see that adding the result (3.6) to the free result (3.2) is 
compatible to the order calculated with the form 



*V(tOord = - v^) Vn m , (3.7) 

where = A(l + 0(e)) and fj n is the dimension of the surface operator n (x) = <9i0(O, x). 
The expression (3.6) implies the one loop corrections 

, ~ N + 2 l7l iV + 2 , nn . 

^ 2 =rf - 2+ iV + 8 £ + '-- ' r ln = ±d-± WT ^e + ... , (3.8) 

in accord with long established results. The subscript ord in (3.7) denotes that this func- 
tion is appropriate to the ordinary transition in a statistical mechanical context. It is 
straightforward then to verify that the functional form in (3.7) is consistent with the lim- 
iting behaviours given by (2.10), since (p 2 is here the lowest dimension operator appearing 
in the operator product (fxf) beyond the identity, and also with (2.13) since in the Dirichlet 
case (j) n is the leading boundary operator appearing in the BOE of 4>. 

For the Neumann case, which is appropriate for the so called special transition, we 
may write 

F^v) sv = C+(l + v^)(l - v"*)*-"* , (3.9) 

where 

^= i d _l_ i^±| e+ ... , (3.10) 

is the scale dimension of the the boundary operator 0(x) = 0(0, x). (3.9) is again in accord 
with the usual OPE for 4>(p and as y' — > or v — > 1 it is easy to see, with (3.1), 

where s 2 is given in (2.12). It is crucial that, with the above results for scale dimensions, 
the next to leading corrections vanish, corresponding to the boundary operator 4> n = 0. 

As v — > the particular boundary conditions are irrelevant so we must require 
^>(0) sp = -P^(O)ord) as exemplified in (3.7.9). By considering the terms oc v 71 ^ 2 we easily 
see, following (2.10), that 

V ' SP =-l-^ + 0( £ 2 ), (3.12) 



which is independent of the normalisation of the operator <p 2 . 

We have also calculated the leading e corrections for the two point function of the 
operator 2 , which represents the energy density in a statistical physics context. There 
are two Feynman graphs shown in figs. (la,b), fig. (la) corresponding to the one loop 
correction to G$(x,x'). 
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figure la figure lb 

Details of the calculation are given in appendix A. After removing divergences the 
result is 

Fp(v)W =2A 2 N^-^e(l + v 4 -(l-v 2 ) 2 ln(l-v 2 ) + 2v 4 ^llnv 2 ±v 2 lnv 2 ) . (3.13) 
N + 8 V 1 — v 2 / 

In the Dirichlet case we therefore have to this order, combining with (3.2,5), 

F, 2 (v) md = Cp ((1 - v^) d -^ + ^e{2v 2 + v 4 ^ In, 2 )) , (3.14) 

for (7^2 = 2A 2 N(1 + 0(e)). In the limits v — > 0, 1 this behaves as 

*V» (v) ord ~ C> - 2Cp ej v^ 2 , (3.15a) 

Fp (v) md ~C 03 ^1 + H^±| e _^(i_ . (3.156) 

The first case corresponds to the leading contributions of the identity and 4> 2 in the OPE 
as expected. For v — > 1 the result corresponds to the boundary operator T(x) = Tn(0, x) 
whose scale dimension is d at any conformal fixed point. In the Neumann case the results 
to this order from (3.13) can be simply expressed by 

F^(v) sv = 4C> (l - e^j v"** + F^(v) old . (3.16) 

For v — > 1 the first term is a constant and so corresponds to the contribution of a boundary 
operator 4> 2 with dimension t)^ 2 = r\$i +0(e 2 ) = 2 + 0(e) which is now present in addition 
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to the operator T. For v — > using F^^) ~ C^2 + C^i^^A^v^^ 2 it is easy to see that 
(3.15a,16) are in accord with (3.12). 

4 O(N) Model for Large N 

Besides the e expansion it is also possible to obtain non trivial fixed points in quantum 
field theories by considering expansions in 1/N for theories with A-component fields [17]. 
Such large A methods have the virtue of allowing calculations of critical exponents at 
conformally invariant fixed points for any dimension, at least in the range 2 < d < 4. The 
simplest example is the non linear 0(A) a model for fields (p a G S N_1 . In a lagrangian 
formalism it is classically equivalent to remove the non linear constraint on <p but intro- 
duce an auxiliary field A, without any kinetic term and with an interaction Cj = \X4> 2 , 
whose field equation then enforces the condition cj) 2 = const.. The 1/N expansion for the 
corresponding quantum field theory may naturally be written in terms of propagators for 
(j), A with a single vertex given by Cj. At the scale invariant fixed point Vasil'ev, Pis'mak 
and Khonkonen [4] formulated the a model in terms of a skeleton graph expansion for the 
two point functions for 0, A which can be solved self consistently for the scaling dimen- 
sions 770, ?7a by imposing conformal invariance in a tractable 1/N expansion. Using these 
methods rj^, rj\, which determine all other bulk critical indices, have been determined to 
0(A -3 ), 0(A -2 ) respectively. To leading order 77^ = \d— 1, r\\ = 2 and results are consis- 
tent with those from the e expansion for the renormalisable ^(7(</> 2 ) 2 , as considered in the 
previous section, with the identification 77A = 77^2 . In recent years the O(N) a model has 
been extensively investigated at its conformal fixed point so that the spectrum of operators 
present and their scaling dimensions are well understood [17,14]. 

Writing the two point functions for the basic fields </> Q , A as 

(r(x)^(x')) = S^G^x, x') , (X(x)X(x')) = (X(x))(X(x')) + G x (x, x') , (4.1) 

since we assume manifest O(N) invariance and {(f) a ) = 0. To zeroth order in 1/N [4] the 
basic equations determining these are equivalent to 

(-V 2 + (A(x)))G^(x, x') = 5 d (x - x') , (4.2a) 

J d d x" G+(x, x"f G x (x", x') = ~^S d (x - x') . (4.26) 

Alternatively in the renormalisable (</> 2 ) 2 theory, with the identification A = \g4> 2 , G\ 
represents the leading contribution at large due to summing all bubble diagrams. With 
no boundary present (4.1a,b) are trivial to solve at a fixed point. If, as required by 
conformal invariance, 

A B 

G4x,x') = ^, G x (x,x') = (4.3) 

then (4.2b) requires in general that 
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while (4.2a) gives = \d — 1, (A) = 0, with A as in (3.3), and hence r]\ = 2 for N — > oo. 
The relation between r)\ and in (4.4) is necessary for the integral in (4.2b) to be 
conformally invariant, using the transformation of the measure in (2.1). It is also useful 
to note that for i]^ = \d — 1 B > for 2 < d < 4 and as d — > 4 B ~ 8e 2 /N. In general the 
normalisations of the fields 0, A are arbitrary but with the above choices the scale of the 
vertex function is determined. Without a boundary, when (A) = 0, conformal invariance 
dictates the general form 

(d> a (x 1 )^(x 2 )X(x 3 )) = S a \ Xij = Xi - Xj . (4.5) 

To leading order in 1/N * 

Crhrh\ 1 



AB 2e ■ < 46 > 

The normalisation independent vertex coupling is then C^\j \A 2 B)^ = 0(N~^) which 
implies formally that the 1/N expansion is consistent. 

The aim here is to extend these results to the situation for a plane boundary. The 
essential results in the large N limit have been obtained by Bray and Moore [10] and also 
Ohno and Okabe [11] some time ago but we differ in making essential use of conformal 
invariance which leads to a a more direct and perhaps simpler derivation of the A two 
point function by solving (4.2b) (in both treatments Dirichlet boundary conditions, or 
the ordinary transition, is more straightforward than Neumann boundary conditions, as 
appropriate for the special transition). Following (2.4) if we write 

G*(x, x') = * f+(t) = —L^F+iv) , F^v) = t^Mt) , (4.7) 
then imposing (4.2a) requires F^(0) = A and also, since rj\ = 2, if we take 



(A(x)> = ^, (4.8) 



satisfies the differential equation 



tf 1 + ^p F * + & + 2 -\A X F^ = 0. (4.9) 



This is easily seen to be soluble in general in terms of hypergeometric functions, the 
appropriate solution with the relevant behaviour as £ — > being 

AF(a,b,-Ae;-£) = A{1 - v 2 ) a F{a,\s -b;\e;v 2 ) , a + b=l, ab=-\A x . (4.10) 



* This may be obtained by calculating the lowest order contribution, using the conformal star 
triangle relation, (<b a (x\)(t>P ' (x2)M x s)) = —5 al3 fd d x — -5- ^t£j to=\ , tj but the result is also 

J (\x\—x)' i \xi—x)' l ) a -l^ (2:3 — X)* 

in accord with a more systematic analysis based on requiring cancellation of shadow fields in the OPE [14]. 
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For the particular applications of interest here we assume that the leading term for in 
the large N limit is given by the iV — > oo results of the previous section. It is important 
to note that the first order in e results for the scaling dimensions are exact in this limit t. 
Hence in the Dirichlet case we obtain 

Mflord = A(£(l + , ^ Word = A(1- V 2 ) ^ . (4.11) 

This solution is in exact agreement with the general form (4.10) for a = \d — 1, b = \e. 
Hence in (4.8) in this case 

A A , ord = -(4-d)(d-2). (4.12) 

The consistency of this result may be verified by considering the OPE for 00 where from 
(4.5) and (4.1,3) in the limit s = x — x' — > 

f(¥(^V^i a ^^rtM, C U X = ^±. (4.13) 

Applying this to the two point function (4.1) with the general form (4.7) then we must 
require as v — > F^(v) ~ A+C^ x A\v 2 . Using the explicit expression (4.11) in conjunction 
with (4.6) is easily seen to recover (4.12). 

Given an explicit formula for then G\ is determined by solving (4.2b). Given that 
rj\ = 2 to leading order in 1/N then conformal invariance dictates the general form 

G ^ x ') = J^h(0 = ^F x (v), F x (v)=eh(0- (4-14) 

We first discuss the essential mathematical problem of inverting an operator rep- 
resented by a kernel G(x,x') which transforms as a scalar of dimension a at x,x' on 
= {(y, x); y > 0}. Representing the inverse by a kernel H{x, x') we require 

/ d d xG(x 1 ,x)H(x,x 2 ) = 5 d (x! - x 2 ) , x u x 2 eR%. (4.15) 

Under a conformal transformation as in (2.1) 5 d (xi — x 2 ) — > 3 (a;i) d 5 d (a:i — X2). Hence 
the integral (4.15) is compatible with restricted conformal transformations preserving the 
boundary if we take 

G ^ = J^9ii) , H(x, x') = j^jfiszzKt) . (4.16) 

Hence we may expect to reduce (4.15) to an equation for the single variable functions g,h. 
To achieve this it is convenient to define a transform, g — > g, by integrating G over planes 
parallel to the boundary* 

Jd d -^G(xy) = J ^ ¥ ^g( P ), P=^tA A=I(d-l), (4.17) 



"I" The two loop results [12,19,20,21] for the surface exponents r), r) n vanish as iV — > 00. 
* This is analogous to the radon transform [22]. 
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where 



g(p) = f(A) Jo duuA_1 0(« + p)- (4.18) 

This transform may be inverted, g — > g, by 

</(0 = n x T{ _ x) J o d P P 9(P + , (4-19) 

where the integral of is singular for d of interest here but may be defined by analytic 

continuation in A from Re(X) < 0. The inversion formula may be verified by using 

Jduip-u)^'^- 1 = A) p£ +A_1 ~T(-\)T(\)6(p) as//->-A. (4.20) 



For <i = 3=^A = lwe use 

-A-l 



P+ <)'{/» for A • 1 . (4.21) 



r(-A) 

to reduce (4.19) to the simple form 

9(0 = --9(0, (4-22) 

which is easy to check directly from (4.18). 

With the definition (4.17) for g, and correspondingly for h, it is easy to see, by 
integrating over x l5 that (4.15) reduces to 

Jp^g(pi)h(p 2 )=4y 1 6(y 1 -y 2 ), Pl = . (4.23) 
If y = e 2d , yi = e 29i this may be simplified to 

/oo 

d6>#(sinh 2 (6> - 6i)) h(smh 2 (6 - 6 2 ))= S(9 1 - 6 2 ) ■ (4.24) 
-oo 

This is then straightforward to solve by Fourier transforms 

~ g (k) = J d6 e lkd g(smh 2 6) , (4.25) 

which finally gives 

l(k) h(k) = 1 . (4.26) 

Hence h may be determined and by taking the inverse Fourier transform and applying 
(4.19) it is possible to find h. 
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For the problem of interest here then, given the result (4.11) for G^, we should take 
a = d - 2 and, with /^(Olrd = ^(0, 

^)=[^+l)]" d+2 - (4-27) 

For simplicity we describe here the application of the above method for the physically 
interesting case of d = 3 when the integrals are also relatively standard. General d is 
discussed in appendix B. For d = 3 we easily obtain 

g(p) =irln^—, g(smh 2 6) = -2tt In | tanh#| , (4.28) 



and therefore 



2tt 2 

g (k) = — tanh iyrA; . (4.29) 

K 



Using (4.26) to give h{k) we then get 

?iW = -iM^W " K) = " ' <430) 

with h — > /i given by (4.22). Hence for <i = 3, when A = 1/4-k, the solution of (4.2b) gives 
in (4.14) 

A(f)o,d = ^ F? ^| F , M- = ^(i-»'). (4.31) 

For general d the results of appendix B lead to 

r(d)r(d - 2) d 1^,00 ^-3(1 + 20 

A(0ard - B r(M _ 4) £ F(d - 2, d, 2d - 4, --) - 2S ^ + , 
16 eT(d-2) 



(4.32) 



iv r(2-|d)r(|d-i)3' 



where Q^_ 3 is an associated Legendre function and, given A in (3.2), B is in accord with 
(4.4) for 7]^=\d- 1. For d = 3, using F(l, 3; 2; z) = |((1 - z)" 2 + (1 - z)" 1 ), this is 
easily seen to be in agreement with (4.31). It is also of interest to consider the limit as 
d — > 4 when 

r d + 2 F(d - 2, d; 2d - 4; - i) = " 2 ' ^ M " 4 ? "|) 

1 + e In £ - 6, £ ^-hfcr^ 7^ (4-33) 



(1+0 2 V ' ^n(n + 2)(n + 3) £ 

(1 - v 2 ) 2 (l - eln(l - t> 2 )) + e^ 2 + z; 4 ^^ lnv 

[ 1 — v 
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With B ~ 8e 2 /N this agrees exactly with the results in (3.4) and (3.14) in the N — > oo 
limit for A -> ±#*0 2 ~ 87r 2 ^ 2 /AT . 

Manifestly from (4.32) /A(£)ord may be expanded in powers £,~ d ~ n or (2£ + l)~ d_2n for 
n = 0, 1, . . . which, from later results, represents the contributions of boundary operators 
with dimensions d + 2n in the BOE. For £ — > oo /A(£)ord oc as expected for the leading 
operator appearing in the BOE for X(x) being Tn(0, x) = T(x). An alternative expansion 
of f\(€) or d valid for small £ also verifies the consistency of the the leading 1/N expression 
for the A two point function with the OPE for \(x)X(x'). From (B.8) in the limit £ — > 
we get 



/a (£)ord = 5 (1 (1 - (d - 2) (d - 3)C) + J (d - 1) (d - 2) (d - 3) (d - 4) In ± + . . . ) . (4.34) 

The leading operator, apart from the identity, appearing in the OPE for AA is A itself and 
the relevant coefficient is determined from the three point function 

(X(x 1 )X(x 2 )X(x 3 )) = - ° 2 XX " lvx ■ (4.35) 

To leading order in 1/AT [14], with our normalisations, 

C XXX = B"^ (4.36) 

and, with r]\ = 2, we can write 

W')~^ + ^A M , C»» = %i, (4.37) 



which implies F\(v) ~ B + C\\ x A\v 2 , using (4.8) for (A). It is easy to see from this that 
(4.34) is compatible with the result (4.12) in this case. Note that C\\x = for d = 3 which 
is reflected by the absence of at; 2 term in (4.31). 

In the Neumann case if we take the N — > oo limit in (3.9) we get 

U(0 SP = A 1 + 2 l d i , FM sp = A(l + v 2 )(l - v 2 )-^ , (4.38) 



where the leading surface operator has dimension d — 3 which agrees with the N — > oo 
limit of (3.10). The expression given in (4.38) corresponds to the general solution (4.10) if 
a = — |e, 6 = 1 + \e so that in this case 



A A , sp = (4-d)(6-d)= £ (2 + £ ). 
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(4.39) 



From (4.12,39) A\ jSp /A\ jOTd = — (2 + e)/(2 — e) which is in accord with the large N limit of 
(3.12). The solution of (4.2b) to find G\ is more involved and is undertaken in appendix 
B. The conformal invariant function f\ in (4.14) is given by 



fx(Os P = B 



1 6-d T(d)T(d-2) 
3 d-2 T(2d- 5) 



- ^ (1+a] b^) 3jP2( ^ + ^^"^^ ;rf_ ^5 ;_ 4^(i+o ) ( 44 °) 

I Trr(Irf-l) 2 8 , 3 7 t 1 , 

+ r(rf-3)r(irf-f)r(|-irf) (1 + 20 2 u ' ' 5 2a '(i + 20 2j 

The two terms correspond in the BOE to two classes of boundary operators with dimensions 
d + 2n and 2 + 2n, n = 0, 1, . . . respectively. Corresponding to (4.34) we also find for the 
singular behaviour as £ — > 

= B (i(i + („ - 3 )( 6 - + '"-^-g^") 4 """ 6 ' 2 4 +•••)• (441) 

in which the first two terms may also be easily seen to be compatible with the OPE (4.37) 
together with (4.39). 

In both (4.34) and (4.41) there appear ln£ terms which are naively unexpected on the 
basis of the OPE. However these may be understood by considering the role of an operator, 
denoted as A 2 , which has dimension r\\i = 4 + 0(iV _1 ). This may be defined in terms of 
the OPE for AA where we take C\\ = 1 and we assume that, with this normalisation, 
A\2 = Aa 2 (1 + 0(-/V -1 )). According to (2.10) we then expect from the OPE a contribution 
to the A two point function given by 

/aa(£)as ~ A\ (l + \ (2t7a,i - ifcv) ^ In ^ + . . . ) , (4.42) 

where t/a,i, ??a 2 ,i are the coefficients of the 1/N terms in the large iV expansion of the 
scale dimensions of A, A 2 . The first term on the r.h.s. of (4.42) clearly represents the 
disconnected pieces (A)(A) in (4.1), which are 0(1) as iV — > oo, while the second 0(iV -1 ) 
term corresponds exactly to the ln£ terms in (4.34,41), with (4.12,39), if 

(^-*-.4=fclNh (4 ' 43) 

which agrees with direct calculations.* 



* rjx t i and r/ X 2 x are given by (5.29) and (5.30) in the fifth paper listed in ref. [18] where 771(5) = 77^1 
with B = 4d(d - 2)77^1 /iV. 
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5 Two Point Functions for the Conserved Current and Energy 
Momentum Tensor in the Large N Limit 

In the O(N) model described in the previous section other operators may be defined 
in terms of the basic fields (fi a through the OPE. Here we focus on the conserved current 
J^P = — J^ a , whose charges generate the O(N) symmetry, and the energy momentum 
tensor T^„ in this model. These have scale dimension d—1 and d respectively. The former 
may be defined through the OPE 

~ ^A- S(iJ « V) , ( 5.i) 

where the coefficient is determined by O(N) Ward identities [14] with Cj setting the scale 
of the two point function for J^f . This may be written as 

(jf(x)jf(x')) = (S^S^ s -S aS S^)Gj^(x,x'), (5.2) 
where, as in (2.16), conformal invariance dictates 

Gj^(x,x') = —±—(l^(s)Cj{v) + X li X' v D J (v)) , Cj = CM . (5.3) 
\ s ) 

In the N — > oo limit rj^ = — 1 and A is given by (3.3). The leading contributions to 
the two point function are then simply 

Gj^(x, x') = (-^j 1 ) \ {d^G^x, x'^'uG^x, x') - d^G^x, x') G^x, x')t' v ) . (5.4) 
Substituting the conformally invariant form for G^ in (3.1) then gives 

CM = (^p)^ ((d ~ 2)FM 2 " (1 " v 2 )vFMFAv)) , 

DM = { F Mv^((l ~ v 2 )vF,'(v)) -Al- v 2 )F^{vf) , 

which leads to Cj = 2/(d— 2)S%. Using (4.9) we may verify that the general results (5.5) 
satisfy the conservation condition (2.18). For the Dirichlet case* using (4.11) we have 

Cj(v) ord = Cj(l + v 2 )(l - v 2 ) d ~ 2 , CMord + DMord = Cj(l - v 2 )^ 1 . (5.6) 



The calculation of two point functions involving the energy momentum tensor 7) 
is more involved. In a non trivial conformal field theory the natural definition of 7) 



* For a free complex scalar field which has a conserved current = d^,4> then in (2.16) using 
(3.2) gives the corresponding results, if Cy = T{±d) 2 / {2-K d {d - 2)), C(y) = C v (l ± v d - 2 (l + v 2 ) + v 2d - 2 ) 
and C{v) + D(v) = C v {l±(d- l)v d - 2 (l - v 2 ) - v 2d - 2 ). 
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is through its appearance in OPEs. In the simplest case we may extend (4.13) to give 
generally 

<f> a (x)<f> a (x') ~ -P^- + ^ C^°(s, d x .)\(x') 

(5.7) 

c% NC^ 1 



d-1 C T Sd ( s 2)^-Jd+i s A* s «' T ^(^) + --- 



where is the coefficient of the bulk two point function for <fi, in (3.1) = F^O) while 
Ct is given by the scale of the two point function for the energy momentum tensor, as in 
(2.32). The coefficient of the energy momentum tensor term in OPEs is determined through 
Ward identities [3,14]. In the term containing A in (5.7) C Vx,0 (s, d) = l+0(s) is introduced 
to include all derivatives of A in the above OPE. It is determined by the requirement of 
reproducing the three point function (4.5) where in general = C\. It is necessary 
to consider such derivative operators, unlike in (5.1), since T^, which has dimension d, is 
not the lowest dimension operator appearing in this OPE. Explicitly 

1 r 1 

C a > b ( s ,d) = — / daa^-^l-a) 8 - -1 

S(a + ,a_) J 

(5.8) 

xVi \-\s 2 a{l - a)d 2 ] m e as - d , 

£- m! {a + l-\d) m 

for a± = |(a ± b) with the Pochammer symbol (a) m = F(a + m)/T(a) and we assume 
[s, d] = to move all derivatives to the right. For our purposes we need to expand this to 
0(* 2 ) 

C a >°(s, d) ~ 1 + \s-d + ((« + 2)s,s v d,d v - - - \ a _ s 2 d 2 ) . (5.9) 

Initially we focus on determining the two point function (T fJ/l/ (x)X(x / )) which as a 
consequence of (2.19,22) has a unique functional form in the conformally invariant limit. 
In order to use the OPE (5.7) to find this we need to determine an expression for the {</><j>\) 
three point function. To leading order in 1/N the connected three point function in the 
O(N) a-model with a plane boundary discussed previously is simply given by 

( ( p a (x 1 ) ( p f3 (x 2 )X(x 3 )) conn = -5 a f [ d d rG 4> (x 1 ,r)G 4> (x 2 ,r)G x (r,x 3 ), (5.10) 

where G^, G\ are given by the results of the section 4 for the ordinary and special critical 
points. The application of the OPE to this is based on the equation (verified in appendix 
C) 

(5.11) 



2m 2rj 2 ~ K"L*,"X2J 2 V 



Y{l d - Vl )Y{\d- m )Y{ V -\d) 
V = Vi+V2, p = ni — r 

r(7/i)r(7/2)r(d-7/) 
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Neglecting the 5 function term, this result is the essential condition which determines the 
derivative operator C a ' b (s,d) in the OPE for scalar fields in conformal theories. For our 
purpose we take r\\ = rj 2 = \d — 1 and generalise this to write 

G ( j } (x 1 ,r)G (j) (x2, r) ~ C d - 2 '°(x 12 , d X2 )G^x 2 , r) 2 + pA 2 (x 2 2 )^C 2 '°(x 12 , d X2 )8 d (x 2 - r) 

- ^(x 12 ) fl (x 12 ) 1 yt^(x 2 , r) + . . . , (5.12) 

where other terms are less singular and irrelevant here (for free fields without boundary 

(5.11) of course shows that and other higher order contributions in the expansion 

(5.12) vanish). Using (5.12) in (5.10) shows that it is compatible with the OPE (5.7), 
taking r)<f, = \d — 1, r\\ = 2 and = A, since the G ( j ) (x 2 , r) 2 term is absent due to (4.2b). 
In this case it is necessary that pA 2 = — A/ (2s) = —C^ x to leading order in 1/N, which 
is in accord with (4.6). If we use the result for Ct when N — > oo, which is the same as for 
free scalar fields [3,14], 

Nd 

(d^Tjs 2 

we may also obtain 

(T^(x)X(x')) = -N [ d d rt^(x,r)G x (r,x'). (5.14) 
Using (5.9,11) and the basic Green function equation (4.2a), with (4.8), gives 

t/j,v(x, X ) t)j,v(%i % ) ~~,$ii,v TZ \oG (f)(x , X ) , 



C T = - J iNo2 , (5.13) 



d ^ (2y) 2 

d 1 

fiuGfj) + _ yVftu {G 2 ) , = d^d v - 



(5.15) 



Using the form (4.7) for G§ becomes 

(2v') 2 / 1 \ 

t^(x,x') = y -j^-[x^X u - ^)f(v) , 

m = - + 1) {(* - 2)f*± (e±F,) - d e(±F 4 



(5.16) 



where is defined as in (2.14). Under conformal transformations t^ v (x,x') transforms 
as a scalar with scale dimension d — 2 at x', and also as a symmetric tensor of dimension 
d at x, so that the integration in (5.14) preserves conformal invariance. From (5.15), or 
from (5.16) and using (4.9), we may also find 

(^d-2)n v + yd v y jG 4> (x,x') 2 , (5.17) 

which in turn gives 

2A\ 

d^t^(x,x) =n v j—^G 4> (x,x') 2 , (5.18) 
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where it is important to note that since f(v) oc v 4 as v — > and hence t fIU (x, x') = 0(s 4_2d ) 
there are no 5 function contributions. Since is easily found from (5.15) or (5.16) and 
using (5.18), and applying the basic equation (4.2b), we therefore can derive 

2 A \ A A 

(T^(x)X(x')) = -J^S d (s) , d,(T,„{x)\{x')) =n,j^ S d (s) . (5.19) 

This is in exact agreement with (2.23) for O — > A, rj — > 2 if we take C = — 2/d. As a 
consequence of the general conformal invariance results given in (2.19,22,26) the integral 
(5.14) is completely determined therefore by 

(T^x)X(x')) = ~j^g- d ^^{X,X V - l f, v )v d . (5.20) 

The same procedures may also be applied to find the energy momentum tensor two 
point function if we start from ((f) a (xi)(j) a (x2)(f) l3 (x 3 )(f) l3 (x4)) and consider the OPEs for 
£12, £34 ~ 0. The relevant contributions for large N are 

N(G t j,(xi, x 3 )G < j,(x 2 , x 4 ) + G^(xi, x 4 )G (j> (x2, x 3 )) 

+ N 2 [ d d r [ dV G^xu r)G4x 2 , r)G (x 3 , r)G (x 4 , r) G\(r, r') . (5 " 21) 

For X12 = £34 = this vanishes due to (4.2b). This is crucial in cancelling the 0(1) 
contribution as x\ 2 ~ £34 — > 0, which would correspond to a 'shadow operator' of dimension 
d — 2 in the OPE (for free fields this is represented by the operator (p 2 ) so that the leading 
behaviour is oc {x 2 2 )^ S i ( x 34)^ s which is appropriate for A, with dimension 2, being the 
lowest dimension operator apart from the identity. Following a similar discussion to the 
above we may now obtain the large N expression for (T fJ , u (x)T (T p(x')) as 

G^p jl>ap (x, x ) — N— {5^ v (2y'Y^ ap ^ X ' X ^ ^ ap (2y) 2 ^ pu ^ X ' X 



2 4> 2 r r (5.22) 

N^S ltv S ap7T ^G <l> (x, x') 2 + N 2 / d d r / dV t^(x, r)t ap (x', r') G x (r, r') , 

""' ) JR d , JR d , 



d 2 "^ P (4yy 

where, with the operator Vp V defined in (5.15), the 'free field' part G? is given by 
G^vap = N^Gff, T> p U T>' ' a P G '0 + T>p V G§ T> ' ap G$ 

d —— (v^ (G^V'vpG^) + V' ap (G^D^G^j ) (523) 



+ 



2(d- 1) 
d 2 



8(d- 1) 



Y)2 v »v v 'vp( g 4?)} ■ 



It is easy to verify that the trace on \xv and op in (5.22) vanishes and hence the terms 
involving A\ explicitly in (5.22) and also contained in t^, t ap can be dropped. Hence we 
can write 

) , (5-24) 
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where 

G\ l/ap (x,x') = N 2 [ d d r [ d d r'i^(x,r)L P (x\r')G x (r,r'). (5.25) 
Writing r M = (z, r) and using the result (5.14,20), 

(5.26) 

for s 2 = (a; — r) 2 , -D 2 = s 2 /(s 2 + 4yz), s /2 , -D' 2 similarly defined with a; — > x', and 
X M (a:,r), X'^a:', r) vectors formed as in (2.14). 

To verify the conservation equations we may use the definition (5.23) along with 
(4.2a,8) to find 

d»G^„ ap (x, x') = ^ ^ ((d - 2)n v + yd v y ap (x', x) , (5.27) 

and using (5.17) in (5.25) it is easy to verify that this is cancelled by d lJL G x lJLV(Tp . By virtue 
of conformal invariance G^ uap (x,x') may be expressed in the form (2.27) but (5.27) now 
gives 



(v± - d)^{v) = jj^aHv) + (d+ ^ d 1 2 V (.) + 2NA x 1 -f(v) , 



(5.28) 



(y± - d)af(v) = 2(d - l) 7 /(«) + 2NA x ^ r (vf'(v) - 3-^/(1;)) • 
By explicit calculation we find using the results in (4.11) or (4.38) 

/Word = 2A 2 ^£ z; 4 (l - v 2 ) d ~ 4 , (5.29a) 
/Wsp = |4l ^ - ^) d " 6 {(^ " 3)(d - 4)(1 + z; 2 ) 2 - d(l - t; 2 ) 2 } . (5.296) 
For the ordinary transition then (5.23), in conjunction with (2.27,38), gives 
a / Word = _ ^ + _ „y-V + 8^(1 - ^"V} , 

y ( „w = -^(^^{^(i-.r-' + ^d-r-v}, (5.30) 

It is easy to check that (5.29a) and (5.30) satisfy (5.28). The evaluation of the integral in 
(5.26) is discussed in appendix D. The results can be similarly expressed in the basis given 
by (2.27). 
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From (D.33.34) we obtain 



(1 - v*) d -\F 2 (1, d - 1, id - 2; d - |, id; -^^) 
d (l-^ 2 ) 2 - 



+ ^( 1 + ^ 1 ^)' 



x 

d 



3 F 2 (1, d - 1, |d - 2; d - |, id - 1; -^^) ( 5 ' 31 ) 



«(l) or d^- 2 (l-^ 4 ), 



4(d- l) 2 



X / \ -W d — 2 2/ 2\ d - 2 



x 3 F 2 (1, d - 1, id - 2; d - i, id - 1; - ^jf^ ) 



+ 4(d-l)'(d + l) a ^^ 1 -^ a » 



where 



^ 2 d- 2 r(|d) 2 r(3-id)r(d- |) 
5 d 2 r(d-i)r(id+i) 



«(i) ord = ^- 2 ^ 2 ; ; v : *r;^ 2; • (5.32) 



These contributions vanish if d = 4 when a(l) rd = 2N / ' S% as expected since the integral 
(5.26) contains a factor 4 — d in A\ and in agreement with the conformal invariant theory 
becoming then of course just a free scalar field theory. For the non trivial case when 
2 < d < 4 it is crucial that adding (5.31) to (5.30) cancels the (1 - v 2 ) d ~ 4 , (1 - i> 2 ) d " 2 
terms in a? (v) and also the (1 — tj 2 ) d-3 term in 7^ (u) so that it vanishes as v — > 1. The 
final result for a can be expressed more compactly as 

a(v) old = (1 - , 2 ) d {l + (1, d, id - 1; d - i, id + 1; -^^) } 

+ a(l) ord ^(l + ^^^). (5.33) 

We have verified that expanding this to 0(e) gives the same expressions as the e expansion 
results in I in the large N limit. When d = 2 the first term in (5.33) vanishes leaving only 
the contribution proportional to a(l) ord = 2iV/5' 2 2 . 

In order to consider how this result for a(v) or d behaves in the limit v — > 0, which is 
appropriate for the OPE, we may use the inversion formula for generalised hypergeometric 



25 



functions [28] to write this alternately as 
a{v) ord = -^2(1 -v ) 



'd 

x 



{ 1 + 2(d- 1) (1 XV ^ S - ^ 1 - ^ 2 - 3 - ^ - (T^p) } 

5 d 2 r(d+i) d + 2(i-v*)*^' 2 ' + *' (i- v 2)* h 1 J 

6 Operator Product Expansions 

As described in section 2 there are constraints on the two point functions in the 
presence of a boundary, which in the conformal limit depend on a single variable function 
of £ or v, arising from the OPE. The coefficients appearing in this are a property of the 
conformal theory on R d without any boundary and dependence on boundary conditions 
arises solely through the one point functions of those operators appearing in the OPE. 
Since, as remarked earlier, only scalar operators have non zero one point functions, which 
have the simple form (2.7), we need only consider for our purposes the contributions of 
such operators to the OPE. In this section we determine the functional forms for the two 
point function in the neighbourhood of a boundary arising from all derivative operators 
formed from a given quasi- primary operator. 

We initially consider the simplest case of just scalar operators when the general con- 
formally invariant three point function without any boundary can be written as 

^-^123 

{0 1 (x 1 )0 2 (x 2 )0 3 (x 3 )) BO bouodar, = (xi22) , ( , 1+ , 2 -,3) ( ,2 3) i( ?72+ ,3-, 1 )(,| i) ^,3 + , 1 -, 2 ) • 

(6.1) 

If C3 denotes the normalisation factor for the two point function of the operator 3: so 
that in general {Oi{x)Oj(x')) = 5ijCi/s 2r] \ and C12 3 = C123/C3 we may rewrite the OPE 
(2.9) to include all derivatives of O3 as 

1 (x)0 2 (x') = ^^^-^ C^-(s,d xf )0 3 (x') , s = x-x', rj_= Vl - V2 , (6.2) 

where C a,b (s, d) is the derivative operator defined in (5.8). By virtue of (5.11) the OPE 
(6.2) exactly reproduces the three point function (6.1) for non coincident points. Of course 
the OPE of 0\ and O2 in general contains contributions from infinitely many quasi-primary 
operators but such a summation is left implicit here. 

For subsequent use it is convenient to generalise the expression (5.8) to 



C a /{s,d) = — / daa^-^l-a) — 1 

B{a+, a_J J 



x 



(6.3) 

T A 7 — r-r- \-h 2 a(i - a)d 2 ] m e as - 9 , 
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which reduces to (5.8) when n = 0. The essential result for application of the OPE in the 
presence of a boundary is 



Cn ,b (s,d^)—— = j— — — — — F(a + ,a-;a + l- n- id; -£) , (6.4) 
(2y') a (2y) a +(2y') a - 2 J 

which is easily obtained by application of the explicit form (6.3). 

Applying (6.2) to (Oi(x)02(x')) then gives in (2.8) the equivalent results for the 
contribution of all derivative operators generated from (9 3 

/i2(0o 8 = Ci2 3 A 3 C^ m+V2 - m) F(±( m + 77-), 1(773 - 77-); 773 + 1 - \d- -C) , 
F 12 {v)o z = Ci2 3 A 3 v^-F^rjs + 77-), ±(773 + 77-) + 1 - \d; 773 + 1 - |d; v 2 ) . 

The leading singular piece as £,v — > of course coincides with (2.10). When ?7_ = and 
?7 3 = d — 2 then -^12(^)03 = 1. This is appropriate for free scalar fields with O3 — ■> </> 2 . In 
this case the above results show that the free field expression for (cf) a (x)(f) l3 (x')}, as given 
by (3.1,2,3) for either Dirichlet or Neumann boundary conditions, is reproduced solely by 
the contribution in the OPE of the identity and the operator 4> 2 and its derivatives, taking 
Cf/,^ = l/N and with as in (3.5). 

We may also derive the corresponding formulae for two point functions involving 
the energy momentum tensor T^ v . The OPE is determined by knowledge of the three 
point function and detailed expressions in the relevant cases for conformal field theories 
in any dimension d were obtained by one of us recently. The construction given de- 
pends essentially on the result that for three points x\,X2,x$ it is possible to construct 
vectors Xi, i = 1,2,3 which transform homogeneously under conformal transformations, 
Xj M — > Q(xi)'R,f ia (xi)Xi a , at each Xi. For X\ we have 

2 

x ^--^r~-^r^ A i - 2 „2 ' ( b - b > 

X 12 X 13 x 12 x 13 

while X 2 ,Xz are obtained by cyclic permutation.* A crucial result is that under 'parallel 
transport' by the inversion transformation, as defined in (2.3), 

2 

x 

I/j,a{ x ij)Xj a = ^-Xin , k 7^ i,j . (6-7) 

x jk 



For the three point function of the energy momentum tensor with two scalar fields O 
of dimension 77 we may then write [3] 

it ( \n( \n( w Vd C ( X^X lv 1 \ {X?)* d 

(T^(x l )0{x 2 )0(x 3 )) no boundary = ~~T~^ J J ^ > ( 6 ' 8 ) 



* It is perhaps worth noting that if x\ — > x, X2 — > x' , £3 — > x, for x^ = (— y, x), then X\ — > X/(2yv) 
and X 2 -> -(1 - f 2 )X7(2y'w), with as given in (2.14). 
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A^( Sl d) = ^--J^ + 0(s). (6.10) 



where the normalisation is determined by Ward identities which relate this three point 
function to the two point function for O. The result (6.8) determines the form of the OPE 
of and O which can be written as 

T»„{x)0{x') = ~^j- d J^1 A ^ d *')°(x') . ( 6 - 9 ) 

where it is easy from (6.8) in the limit x\2 to see that 

1 

Using the full result (6.8), with the definition (6.6) for X 1 , we may find after some work, 
in terms of C"' 6 (s, d) given by (6.3), that this may be extended to include derivative terms 
in the form 

A^(s, d) = y 2 d^C^ +2 ^( S , d)-- d 6^C?> d -*( s , d) + terms oc s „ or s v . 

(6.11) 

In the presence of a boundary the OPE (6.9) gives 

(T.Jx)O(x')) = --P^^—L^A u Js,d x A 7 -^, (6.12) 
and applying (6.4) with the explicit expression (6.11) leads to 



-^^(H,-^,-^;^) + termsoc^or^} 



(2y) 



(6.13) 



where in the last line we have completed the s p , s u terms to achieve the desired general 
expression dictated by conformal invariance in accord with (2.19). Clearly this result agrees 
exactly with with the general expression obtained in (2.19,22) with also the result (2.26) 
for the coefficient cto- Of course getting the correct form is a necessary consistency check 
of the above treatment. 

In order to determine the OPE for two energy momentum tensors involving a scalar 
operator O we need the three point function [3] 

no boundary 

1 ( x 2 2 ynf ( X lp X lu 1 \ ( X 2a X 2p 1 \ 

= teF^PP l w(ll2)Co + \~icf- - dH \~icf- - f° p ) Ao 

+ (x? 2 (X lp X 2a I up (x 12 ) + n++v,a <-> p) 



(6.14) 



4 X 2a X 2p 4 X^X lv 4 x \ 
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where Aa^o^o are constants and X Mi , )(Tp is defined previously in (2.28). The three 
terms in (6.14) are separately conformally invariant. However imposing the conservation 
equation for T pv , using (2.29) and also 

d^X 2a = ±-I, a (x 12 ) , =-(d-l)j^, 
di^^jJrjd(x? 2 (X lfl I ua (x 12 ) + X lu I^(x 12 )) + 2 6 ^J[f)) ( 6 - 15 ) 
= - 7 2^ 2 \ ( d ~ 2 ) f (d + l)La(x 12 ) + 2xl 2 X lv X 2a ) , 

\ x 12) X 13 a V J 



gives rise to the conditions 



\-qdCo + A + (cT - 4 - 7]d)B = , 
(1 + \{d -v)(d- l))Ao + (v + 2)(d - 2)£ = . 



(6.16) 



In consequence there remains a single constant parameterising the three point function 
(6.14) but there are no Ward identities which relate this to any two point amplitude in 
this case. 

From the three point function (6.14) we may derive the contribution of the operator 
O to the OPE for two energy momentum tensors which can be written as 

T^(x)T ap (x') = ]_ A^ ap (s, d x ,)0{x') , (6.17) 

where A^ uap (s, d) is determined by requiring it to reproduce the full expression on the r.h.s. 
of (6.14). In the presence of a boundary the operator O then gives rise to an expression 
for the two point function {T pu (x)T cr p(x')) using the result (2.7) for (0(x')). Following a 
similar, albeit more tedious, procedure to that outlined above for (T tlv (x)0(x')) the result 
is compatible with the form (2.26) where the corresponding form for the invariant functions 
A, B, C due to the operator O in the OPE is given by 

+^(4 +2 ^{o>«), 

C(v)o - £&>(Co + Bo A ? iL + Av-^e^HO , 
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We may solve (6.16) by writing 

d-2 

A = C T toj—^v(v + 2 ) , 

Bo = C TT o - \{d + 2 - V )) V , (6.19) 

Co=C TT o(-^ + (d-v) 2 ), 

and then, using various properties of hypergeometric functions and with the results in 
(2.36,37,38), we may obtain, with Ctt° = Ctto I Co-, 

a(v) = C T T°A (d - 2)(d + 1)^(1 - v 2 ) 

{('-^K^^^^W)}. 



x 

l{v)o = C T T A ^^{d + 1K(1 - v 2 ) (6.20) 



X 



2^ 



F(v 2 ) = F(±r],±ri + 1 - ±d;r? + 1 - \d;v , 

Although these expressions are somewhat lengthy they satisfy differential equations equiv- 
alent to the conservation conditions (2.31a,b). 

As a check on (6.20) we may consider the trivial case when O — > 1, the identity 
operator, with ij = 0. In this case in (6.14) C\ = Ct, as defined in (2.32), and A\ = 1 and 
from (6.19) C T = C T Tid(d - 2){d + l)/(d - 1). Using F(t> 2 ) ~ 1 - i?7ln(l-i; 2 ) for r] ~ 
then taking the limit 77 — > in (6.20) and expressing the coefficient in terms of Ct gives 
the expected results for 0:1,71 independent of v. 

If rj = d — 2 then F(v 2 ) = 1. As discussed earlier this case is relevant for the operator 
(ft 2 in free field theory when now Ctt^ = \d(d — 2) 2 A/(d— 1), = ±AN, with A given 
in (3.3). In this case (6.20) becomes 



«(<V= ± N J 2 {d-2){d+l){l-v 2 ) 2 v d - 2 , 

(6.21) 



7(*V = =F 4(d "* 2S 2 (d ~ 2)(d + 1)(1 - t; 4 )^- 2 , 

which is identical to the results calculated in I for free scalar field theory for the terms in 
(T /jll/ (x)T a p(x')) dependent on the boundary conditions. In general however the expression 
obtained from a single operator O in the OPE, as given by (6.20), fails to satisfy the 
required behaviour as y, y' — > 0, in particular boundary conditions such as 7(1) = 0.* 



* In two dimensions the contribution of general scalar operators from (6.20) is zero. Nevertheless 
the two point function (2.34) is reproduced by the OPE T zz (z)T 2Z (z') = e {z - z ") d * 0{x') in terms of the 
quasi-primary operator 0{x) = T zz (z)Tzz(z), which has dimension 4 and Aq = \Ct- 
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In the O(N) model for large N the result for a(v) in (5.34) shows that the scalar 
operators appearing in the OPE of two energy momentum tensors have dimensions 2 + 2n 
and also 2d + 2n for n = 0, 1, . . .. 

7 Boundary Operator Expansion 

In the previous section we described how the contribution of all derivative operators 
formed from a quasi- primary operator and appearing in the OPE for operators 0\ , Oi 
to the two point function (0\(x)02(x')) in the presence of a boundary may be explicitly 
calculated. In this section we show how this can also be achieved for all derivative operators 
formed from a particular boundary operator O appearing in the BOE of 0\ or Oi. We 
assume a basis of boundary operators which have a well defined spin under 0(d — 1) 
rotations and scale dimension under scale transformations. For O a scalar operator of scale 
dimension fj, so that the corresponding derivative operators have dimension fj + n, n = 
1,2,.. ., the two point function on the boundary has the form 

(6(x)(9(x')> = , s = x-x'. (7.1) 

The contribution to the BOE in (2.11) for a scalar operator O of dimension r\ involving 
derivatives of the operator O may then written as 

° {X)= (2^^ (2/2 ^ )a(x) ' (7 ' 2) 

where the differential operator D f '{y 2 V 2 ), V» = <9;, is determined by consistency with 
(2.12). Defining B ° = B o6 /C d this requires 



m=0 

It is easy to see that this is satisfied if we take 



^V^E^^T^i-yvT. (7.4) 



For application to the BOE of two operators at points away from the boundary we 
extend the result (7.3) to 
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Applying this to {0\(x)02{x')) we can extend the leading order result (2.13) to give 

/ 12 (0a = B 0i6 B 02 6 ^F(rj,rj + 1- \d- 2fj + 2 - d; ~) , (7.6) 

which therefore includes all derivatives of O in the BOE of 0\ and 2 - 

As a special case we may consider fj = \d — 1 when F(fj, A; 2A; z) — > \ ((1 — z) - '' + l) 
as A — > 0. This is relevant for free scalar field theory with Neumann boundary conditions 
for the surface operator 0(x) = 0(0, x) when B^ = 1 and (7^ = 2A so that (7.6) gives 

which is the exact result in this case. Alternatively if 77 = |<i then the hypergeo metric 
function simplifies to F(rj, 1; 2; 2) = ((1 — z) 1-11 — 1) /{fj — l)z. This is appropriate for free 
scalar field theory with Dirichlet boundary conditions for the surface operator n (x) = 
#10(0, x) when now B^ n = \ and Cg = 2(d - 2) A so that (7.6) becomes 

M£W = - (1 + ^d-i) ' ^ 

which is again of course the exact result. 

Another case when simple formulae are obtained is when f) = d — 2. This is relevant 
for the boundary operator 2 (x) = 0(0, x) 2 in free field theory with Neumann boundary 
conditions. Applying this to the two point function of 2 , with B^k = 1 and C^ 2 = 8 AM 2 , 
we then obtain 

= 8NA " vd ~ 2 > F 4>< V U = + 2NA 2 (1 - v d ~ 2 ) 2 , (7.9) 

where -F^2(i>) sp is given by (3.2,5). The remaining part, after subtraction of the contribu- 
tion arising from 2 , is identical with the result in the Dirichlet case for which the leading 
behaviour as v — > 1 is produced by an operator of dimension d. 

We may also extend this treatment to the BOE for operators with spin although 
we confine our attention here to the energy momentum tensor T^ v . Following (2.14) for 
x^ = (y, x) and x' defining a point on the boundary we may define a vector under conformal 
transformations by 

= t|s m -rift = -I^ui^riu, X 2 = l = (y, x - x') . (7.10) 

s 

Using 
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it is easy that T Mi , can only have a non zero two point function, invariant under conformal 
transformations and satisfying the conservation equation, with a scalar surface operator 
provided that this operator has dimension d. In any conformal theory there is such a surface 
operator T(x) = Tn(0,x), since Tn is non singular as the boundary is approached, and 
it is natural to assume that this is unique. The surface two point function for T has the 
form (7.1) with a coefficient Cf. From the general result (2.27) or specifically from (2.40) 
we may then show that 

Cf = ^C(l) + (A(l) + 4B(1)) = a(l) , (7.12) 

with a(v) defined in (2.37,38). The two point function of T^ v and T can then be written 

(T,„(x)f(x')) = -^-Cf ± (X»X V - 1^) , (7.13) 

where we have taken Bt t = 1. 

From (7.13) we may deduce the form of the BOE for T Mi , involving T, 

T^(x) = D^(y,V)T(x), (7.14) 

where we must require 

D *»(v> ^ = j^M*^ - ■ (7 - i5) 

Writing 

D^(y, V) = D^(y, d)D d {y 2 V 2 ) , (7.16) 
with D d (y 2 V 2 ) defined by (7.4), then we may take* 



1 d 

^ = (d-l)(d+l) { y2dl3 l ~ Sijy dy~ ^ Sij ) ' 

Ai(y, V) = + l + y—)d l , Dud,, a) = -Ai(y, 5) 



(7.17) 



With these results we may calculate the two point function for the energy momentum 
tensor T^ v and a scalar operator O of dimension rj by using the BOE involving T which 



* Since (y 2 d 2 + (d + 2)yn-d)l/s 2d = there is some ambiguity in the definition of D^ u (y,d) but 
this does not affect the final results in the BOE. 
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gives 

(T^(x)0(x')) = B f C f (2y') d ^D pu (y^)D d (y' 2 V' 2 ) ' 



s 2d 



Bo*C t W)*-*b^(y,d)j^F(t) 

T /y (V)"" 1 \( v v 1 



- (n M n, - + 1) _ + d (e + 1) - - d) f(o j 

for F(0 = ^F(d,|d+l;d+2;-l). 

It is easy to see that this is in agreement with the general form given by (2.19,22) and 
comparing with the result (2.26) shows that we must have 



B of = B T C f = -V^, (7-19) 



S d 

which was first obtained by Cardy [7] for consistency of the OPE and BOE. 

We also consider the contribution of the boundary operator T to the two point function 
of the energy momentum tensor where using (7.14) we obtain 

(T^{x)T ap {x')) = CfD^y, V)D ap (y', V')-^ 

d . . , * , 1 , ^ (7-20) 



d 



-C t D ap (y\ d>) (— (X»X V - 1 ^) [ft + 1)] 



using the result already obtained in (7.19). The evaluation of (7.20) is straightforward 
albeit tedious. Some details are given in appendix E. The final expression is of the necessary 
form shown in (2.27) where the coefficients are given by 



C(v) f = rr—^n , ~, C f \v d ~ 2 {l - v 2 ) 2 , 
d 



d 

7 " (d-l) 2 (d+l)^ T2 ' 
2B{v) t + C(v) f = jj^C t \v d ~ 2 {l - v 4 ) , (7.21) 

A(v) f + 4B(v) f = {d _ 1 y {d+1) Cf i y~ 2 ((d + 2)(1 + v*) + 2dv 2 ) . 

These results satisfy the conservation equations (2.31a,b). Further for d = 2 they are exact 
since clearly only scalar operators on the boundary need be considered and in this case 
(7.21) coincides with (2.33) if C T = Cf- From (2.37,38) we may also write 

a(v) f = -f^—rCf ±v d - 2 (d(l + v 4 ) + 2(d - 2)v 2 ) , a(l) = C f . (7.22) 
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For application to a non trivial conformal field theory we may consider first results to 
first order in e at the non Gaussian fixed point in scalar field theory. In ref. [23] and in I 
Cf was calculated to this order for Neumann/Dirichlet boundary conditions giving 

N , 5JV + 2 

The coefficient appearing in (3.15b) determines {B,^ T oxs \) 2 Cf ord and hence using (7.19) 
with O — > 4> 2 , together with (3.8) for 77^2, we find 

(V,ord) 2 _ N 



Cf, = —. j 2± — £ . (7.23) 



Crh2 2 



(l + 0(e 2 )). (7.24) 



It is also of interest to compare the expressions obtained in (7.21,22) with the results of 
the 0(e) calculations for (T^(x)T ap (x')) in I. In the Neumann case for v — > 1 



a(v) sp -a(v) f (1-zrf-^ 



/V -, _ iV+2 

7Wsp-7Mf~ " oaf (1-^) ^ , (7-25) 



AT JV+2 



^d 



in terms of the functions defined in (2.36,37,38). These results correspond to the contri- 
butions expected to arise from a boundary operator which is a symmetric traceless tensor 
7lj of dimension f)j- = d— + 0(e 2 ). For such an operator we may define a two point 

function with the energy momentum tensor 

(T^xfaW) = B Tt {2y) ^ (l^aiS) + ^- [ (x,X v - ^<W)%) , (7.26) 

using X which is defined in (2.28). This satisfies the required conformal transformation 
properties for T Mi , and obeys the necessary conservation equation for arbitrary fjf. In 
the Neumann case the boundary operator %j{x) relevant for the limiting behaviour in 
(7.24) appears as the leading term in the BOE of the operator formed from the traceless 
part of Tij(x). The above analysis shows that the coefficient, oc (2y)~ d+fl T , is singular 
for y — > unlike the situation for Tu(x) = —Ta(x) which tends smoothly to T(x). %j 
cannot contribute to the BOE of scalar operators so its role is not apparent in previous 
discussions. In the Dirichlet case, corresponding to (7.24), we have 



a(v)ord 


— a(v)f ~ 




7(^)ord 


-7(v)f ~ 


AT JV + 2 _ 


e(^)ord 


-e(v) f ~ 


AT IV + 2_ 

^d 



(7.27) 
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which represents the contribution for an operator %j with dimension d + 2 — ^§£- This 
is as expected since such an operator should be constructed from the field <fi in terms of 
expressions of the form did\§djd\§. 

In the critical O(N) model in the iV — > oo limit we may also use the results obtained 
in sections 4 and 5 to verify the consequences of the BOE. By considering the limit £ — > oo 
or v — > 1 of the A two point function we may identify the contribution of the boundary 
operator T of dimension d. In the ordinary case, from the 0(£~ d ) term in (4.32) we may 
identify 

[B x ordj C f ord = B v ^ 2d _^ ■ (7-28) 
Using now (7.19), with 77 = 2 and A\ j0r< i given by (4.12), we then find 

~ _2N r(2rf-3)r(3-irf)r(irf) 3 

ford " S% T{d)T{d-iy ' 1 } 

Since a(l) = (7j, this is identical with the result obtained by direct calculation in (5.32). 
The results for the energy momentum tensor two point function given in (5.30,31) and 
(5.33) exhibit explicitly the full contributions of the boundary operator T, as given in 
(7.21,22), while the remaining parts involving 3F2 functions arise from non scalar boundary 
operators with dimension 2d — 2 + 2n. For the special case (4.40) gives 

(Bx s P ) C fjSp - B- d _ 2 r(M _ 5) , (7.30) 

which then implies 

c TsP - -gj 6(6 d) r(d)r(d _i )2 • ( 7 - 31 ) 

It is easy to check that (7.29,31) are in accord with the e expansion results in (7.23). 
8 Conclusion 

Although the critical behaviour of statistical systems with a boundary is relatively 
unexplored experimental investigation is feasible. In this paper we have discussed theoret- 
ically the form of the functional behaviour of two point functions at a conformal invariant 
crtical point in dimensions d > 2. In particular the O(N) model, for iV — > 00, provides 
a tractable non trivial example which may be analysed for 2 < d < 4. The results of cal- 
culations in this model, which have been here extended to include correlation functions 
involving the energy momentum tensor, are complementary to the £ expansion and provide 
a limiting case for other approximations. 

We may perhaps note that the functions of the invariants v or £ are initially defined 
on the physical region 0<t><lor0<£<oo but may be analytically continued to the 
whole complex plane. The singularities as v — > and v — > 1 are well understood in terms 
of the OPE and BOE. It would be interesting to understand more directly the form of the 
singular behaviour asv— >ooor£— >— 1. In this context we may note that our results 
have simple transformation properties under v — > v -1 or £ ^ — 1 — £ (which is equivalent 
to taking y — > — y or y' — > —y') which may merit further investigation. 
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Appendix A 



Here we present some details of the perturbative calculation of the two point function 
of 4> 2 to first order in e = 4 — d assuming the renormalisable interaction £ int = ^gg{4> 2 ) 2 
with (p a an iV-component scalar field. The coupling is restricted to the non-gaussian fixed 
point g = g* = 487T 2 e/(N + 8) to lowest order. To simplify the calculation we choose 
a perpendicular configuration where = (y, 0) and x'^ = (y',0), with y > y', so that 
v = (y — y')/(y + y'). For the free field case the basic propagator from (3.1,2) in this case 
is just 

To first order in g we may write 

(</> 2 (y, 0)<f> 2 (y', 0)> = 2NG(y, y'f =F |iV(iV + 2)^G(y, y') /« - |iV(iV + J 6 , (A2) 

where J a , 7, are integrals corresponding to fig. (la,b). To this order we may restrict these 
to d = 4 since g* oc e and 7 a then has the form 

U-A'ftJ*,^^)^^. (A3) 

withX 2 =x 2 + (y-^) 2 , X' 2 = x 2 + (y'-z) 2 ,X 2 = x 2 + (y + z) 2 , and X' 2 = x 2 + (y' + zf. 
For finiteness we exclude from the ^-integral an e neighbourhood of the boundary z = 0. 
In the Dirichlet case the result is finite for e — > 0, but in the Neumann case the divergence 
must be removed by the addition of a surface counter term oc (p 2 . This calculation is 
equivalent to finding the one loop correction to G^(x, x') which was calculated previously 
[9,14]. Consequently we easily obtain the result 

„ 2N(N + 2)A 2 e / 4 (1 - v 2 ) . n 2 , , 2 . (l-t> 2 ) 2 \ , . 

<^ 2 («)« = (iV + 8) [ v ln + ln C 1 " « 2 ) ± « 2 ln „2 J • ( A4 ) 

For lb when (1 = 4 we have 



1 1 \ 2 / 1 1 



The integrals involved here are more difficult than the previous case. They contain the 
usual short distance divergence arising from the singular behaviour of G^x^x') 2 when 
d = 4 but this is removed by hand by restricting the integral over the vertex point (z, x) 
to exclude z from e neighbourhoods of y, y'. Evaluating the integral then gives 

7> (v) b = 2N f N + + 2 i A2£ (l + » 4 + v 2 ln (1 - , 2 ) 2 + , 4 In v 2 + ^ In, 4 

4e 2 x (A - 6) 
± v 2 ln (1 - v 2 ) 2 + (1 ± v 2 ) 2 ln • 

(y-y ) J 
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The divergent term as e — > is proportional to G(y, y') 2 In {{y — y') 2 /e 2 ) which of course 
reflects the modification of the singularity as y — > y' from its free value arising from the 
one loop correction to the scaling dimension of 4> 2 . The complete result to this order is 
now 

(v) = 2NA 2 (1 ± v d ~ 2 ) 2 + (v) a + i> (v) b . (A.7) 
Dropping the divergent piece then gives the result quoted in (3.13). 

Appendix B 

In order discuss the evaluation of the integrals in section 4 for general d when g{£) is 
given by (4.27) it is convenient to work backwards and consider an expression for g(k) of 
the form 

, r(q-ifc)r(q+ifc) 

9a > b[k) - T(b-ik)T(b+ik) ' { } 

Having found the corresponding g a ,b(0 it is then trivial to find the associated fo(£), which 
determines the inverse kernel H(x,x'), since it is obviously equal to <7&, a (0- The inverse 
Fourier transform can be found as a sum of residues of poles, at |fc = a + n for 9 > 0, 
giving a hypergeometric function 

j„, l ,(sinh 2 9) = i F jike- m 'g(k) 

= =jz r \ F(2a, a + b - ±; 2a + 2b - 1; — . 

r(6-a)r(6 + a) (4cosh 2 6>) 2a V 2 coshV 

The transform g — > (7 is then straightforward since we can use 



r(A)j„ (i + p + u)p+» (1 + />)"' 

Applying this term by term in (B. 2), with p = 2a + n, we get 
T(2a + A) 



9a,b(0 



4^%^ - a)r(6 + a) (TRF* F ^ + *, a + 6 - j; 2a + 26 - 1; — ) 
r(2f, + A) 1 F(2a + A,a + 6-|;2a + 26-l;-^). (fl.4) 



4 2a - 1 7r A r(6-a)r(6 + a) £ 2a + A 



For application in the Dirichlet case we take 



7r A+1 r(2a)r(irf-a 



»(0 " " IT- r(a+ l )r(a) 2 ^-1A,1(A+1)(0 , (B.5) 
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and hence in this case 



KO = /(a) r(2 ™ a) ^( d ,a;2a;-l), (B.6) 

where /(a) is as in (4.4). Alternatively (B.6) can be written as 
£ d h(£) = f(a) (v 2a F(a, 2a - d; a - d + 1; w 2 ) 



so that h(£) ~ f(oL)^ a d as £ — > 0, assuming a < d. When a = d — 2, as required to 
leading order in 1/N, then 

= f(d - 2) 1 1 (1 - (d - 2)(3 - d)£) + (d - 4) 4 ( In i §F(d, 5 - d; 3; fl + &(£)) } , 



n=0 



(5.8) 

For the Neumann case we consider functions of the form 

Writing (1 + 2£) 2 = 1 + 4£(1 + £) it is easy to see from the above discussion that the 
corresponding transform becomes 

~ m = ^ r£+-fe? ( " 2+ ^ 2) S --Wrf) • " 2 = A2 - 2 («- 1 ) ■ < B10 > 

Defining, as in (4.26), h(k) = l/g(k) then two approaches are possible in order to obtain 
h(£). Firstly finding the inverse Fourier transform may be simplified to the previous 
discussion by removing the pole at k 2 = —An 2 which is equivalent to expressing h in terms 
of a differential equation, 

(" 2 - w) ^ ^ ^ r£yfK hi^-i^ 0) ■ ("•"> 

which can further be rewritten as 



,d 2 . i, d 2 \ £/ , \t, \, i, 1N T(a-l)r(A+l) 

p(i + p)^f + (p + % - ^ 2 J Mp) = -**/(«)(« - id - i) r( L-2)r ( d- a ) 



dp 2 



x (i+p)^ F(A+1>a ~ 1;2a ~ 2; rb ) 



(5.12) 
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By considering the definition of the transform (4.18), or its inverse (4.19), we can then 
derive an equivalent differential equation for h 

= -/(«)(« - ¥ - Df^l. i_ F(d , a _ 1; 2Q _ 2; _ i) , 

which is essentially identical with an equation obtained, by very different means, by Ohno 
and Okabe [11] when a = d — 2 which is appropriate for discussion of the large N limit in 
section 4. The homogeneous equation is easily solved in terms of standard hypergeo metric 
functions. The relevant boundary conditions, which ensure a unique solution, are that as 
£ — > oo, taking p > 0, /i(£) ~ £ _At_A ,£ _d (with no £^ -A behaviour) and as £ — > then 
/i(£) ~ l,£ Q - d (with no terms oc 

However a more direct approach is to take the inverse Fourier transform h — > /i when 
we obtain by contour integration 

t r • u2 m 4A+1 r(a + i)r(a) 2 / 4r(A + l) . 2 

Msmh A) = ^r r(2a)r( i d _ a) ( v r (a -i)r(«-A-f)(^-(A + i) 2 ) G(smh " } 

r(i(A + ^ + i))r(i(A-^ + i)) i |0| 

r(a - 1 - |(A + //))r(a - 1 - |(A - //)) p e 

(B.14) 

with 

G(sinh 2 0) = e- 2(A+1)|e| 4 F 3 (A + 1, B, C+, C_; a - ±, C+ + 1, C_ + 1; e" 4|e| ) , (fl.15) 
where 

S = A + f-a, C± = §(A+l±/x). (B.16) 

The first term on the r.h.s of equation (B.14) represents a particular inhomogeneous so- 
lution to (B.ll) while the second term is a solution to the homogeneous equation. It is 
possible to take the inverse transform h(p) — > /i(£) by noting that with p = sinh 2 9 then 



e 2aW = (Vp + V^) = 4a(1 \ p)a F(a, a + |; 2a + 1; ^ , (5.17) 
so that using (B.3) gives 

i>TA) f(=A) i d " ( " - WT7r F{a - a + h 2o + 1; TT7» 

= ( 1+; ).+x - F (" + A > a +S' 2o + 1 'Y^) (B.18) 
= ^^l^iyf + A + 1), J(a - A) + l; a + 1; . 
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This formula can be used for both terms present in (B.14) if G(sinh 2 6>), as given by 
(B.15), is expanded in a series of exponentials in 9. The part of the final result for h(£) 
corresponding to the e~ 2 ^\ Q \ term in (B.16) is then quite simple, 

4 2C ~ I> + ^)i» 2 T(C+)T(C-) r(A + /i) 

n d T(2a)r(irf -a)Y{a-\- C+)T(a - ± - C_) Tiji + 1) 

, + i 1 (B.19) 

but finding a nice form for the transform of the piece involving the 4-F3 generalised hyper- 
geometric function poses some challenges, a resolution of which is illustrated below. 

After expanding 4 F 3 in (B.15) and using (B.18) for a — > A + 1 + 2n we may then define 

y, 1 (A+1)„(B)„(C+)„(CL)„ (rf) 2 „ 

(4?(l + 0) A + 1 ,t5 1 n!(a-i)„(C + + l)„(C_+l)„(A+l)2„ 1 ' ' 

x 4..(4^(1 (A + 1 + 1 + " ; A + 2 + 2 " ; "WTfl > ' 

so that, in addition to (B.19), the remaining part of h(£) becomes 

4 A+1 r(a + !)i» 2 r(d) i 

7r d r(2a)r(irf-a)r(a-i)r(a-A-f)C + C_ yq) ' y ' ] 

The result (B.20) for F(£) may be rewritten by expanding the hypergeometric function as 

where is given by the finite sum 

_ (-l) n fN\ (A + l) n (£) n (C + ) n (C-) n (cQ 2n (A + 1 + n)y-n(| + n) N . 
N 4n ^ (a-i) n (C + + l) n (C_ + l) n (A + l) 2n (A + 2 + 2n) iV _, 

(S.23) 

By using identities for the T function, such as (d)2 n = 4 n (|rf) n (A + l) n , and (B.16) we can 
write this in the compact form 

hN = <y2 ^\ <y + 2)J N 7-Pe(A + 1,S,C+,C7_,A+ |, |A+ §, -2V; 

a - |, C+ + 1, C_ + 1, \\ + §, §, A + 2 + AT; l) . 
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-n 
iN-n 



All that remains is to find a closed expression for the sum of the terminating series repre- 
sented by the 7-F5 generalized hypergeometric function. For the present this remains elusive 
for general a. However, if we take a = d — 2, as required for section 4, then \i = |(5 — d) 
so that C + = I and C_ = A — |. Consequently the jFq series reduces to a finite 5F4 series 
which is summable by a special case of DougalPs theorem [27] 

5 F 4 (A + 1, \ - A, A - i, \\ + f , - N; d - f , § , ±A + i, A + 2 + N; l) 

= (A + 2) jV (A-|) jV (5.25) 

(d - §)iv(§)iv 

so that in this case 

m = mllr 1 lMX + A - ^ *" - *■ * ■ {BM) 

Hence, with / defined in (4.4), the final result becomes 

HO = f(d - 2)r(|d - 1){ f (6 - ^^^Tj m 



sin7rAr(A — -|) 8 , 3 1 

+ T^T r(d-3) (1 + 20 2 ^' 5 (1 + 20 2 



(5.27) 



The form (B.27), with (B.26), for h(£) is not appropriate for considering the limit 
£ — > 0. However by using the relation* 

3 F 2 (a, 6, p;c,p + 1; -z) 

r(a - P )r(6 - p)r( c ) 1 



r(p + i) ; 



r(a)r(6)r(c-p) ^ 



r(c)r(6-a) p _„ „ , , 1, (5.28) 

- r(6)r(c-a) ^*~ 3F 2 (a,a + l-c,a-p;a + l-6,a + l-p;--) 

~ w1rr°~L P .-V 2 (6,6+l-c,6-p;6+l-a,6+l-p;-l), 
r(a)T(c — b) b — p 2/ 

we may find an alternative expression for F(£), 

m = m T (d-l)T(X-l)T(X-l) £ + i 



r(rf- 4 ) r(A-i)r(A + i) + 

- T{d ~P r — 2 f 1 + (d - 3)(6 - d)£(£ + 1)) 

r(A+i)r(Jd-i) 6-d [4|(i + o] V ' 



(5.29) 



r(A-l)r(|d-3) d-2 2 ^ 2) 



x ( ln4£(£ + 1) 3F2 (A + 1, 4 - §d, id - 1; 3, |d; -4£(£ + 1)) + ^(0) , 



* This does not seem to appear in standard references but related results were found long ago [28] . 
The particular result may be derived by multiplying by (— z) p and differentiating when it reduces to a well 
know result for standard hypergeometric functions, the constant of integration is determined by taking 
z=l. 
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where is given by a power series in £, analogous to h(£) in (B.8). With this form the 
behaviour for £ ~ is manifest. It can also be readily shown that 



1 >, 1 \ 2v^r(A-i) 2£ + l 



(2£ + l) 2 V2 ' ' '(2£+l) 2 ' sin7rAr(A-2) ^ (1 + e) ]A-| 

+ A^F(A-l,l;A+i;-4£(£+l)). 
A ~ 2 

Hence in (B.27) the first terms on the r.h.s.'s of (B.29,30) cancel and 
HO - f(d - 2) * + \.M + (d - 3)(6 - d)£(* + 1)) 



(fl.30) 



+ fid - 2) (rf - 1)( ^f 2 - 4)(rf - 6)2 (^ + 1) m [4,(, + 1)] - . 



(B.31) 



Appendix C 

In order to justify (5.11), and also determine the appropriate form for the derivative 
operator which is exhibited in (5.8) and plays a crucial role in the OPE, we consider first 
the Fourier transform 

J (r 2 + e ~ m W) K ^- l Am 

= ^ (r(J d - rj) [y) ^ d Fl (rj - \d + 1; I^V) 

+ r(r 7 -|d)( / , 2 )^F 1 (id-r 7 + l;i / ,V)), 
where Fi(a;2) = z n /(nl(a) n ) is defined by series expansion. Hence we may write 

d d r e ip ' r 

(x - r) 2 ^(x' - r) 2r >2 



r [daa Vl - 1 (l-a) V2 - 1 e ip < ax+{1 - a)x ' ) ( d a r — r" 

h,m) Jo J 



J d r piyr 

B(vi,V2)Jo K ' J (r 2 + a(l - a)s 2 ) v 

1 7T^ * 



f daa^-^l - a )V2-i e iap.s+i P .x' 
Jo 



B( Vl , V2 )T( V ) J 

x {r(irf - vHW-^Fi {v-id+ 1; Ja(l - a) a V) (^ 2 ) 
+ F(r, - \d) (a(l - a)s 2 ) (|d - 7/ + 1; Ja(l - «)«V) } 
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where s = x — x' and i] = i]i + 7] 2 and C a ' b (s, ip) is defined by (5.8). It is easy to see that 
this is then equivalent to (5.11). 

Appendix D 

For the calculations in section 5 it is necessary to evaluate conformally invariant 
integrals over Mi. Techniques for dealing with conformal integrals on all M d are well 
known [29] . Here we show how the methods of section 4 can be used to calculate the form 
of the integral in (5.26) although a general discussion will be given elsewhere [30]. For 
x = (y, x) and x' = (y', x') the integrals to be considered are of the basic form 

(D.l) 

where restricted conformal invariance guarantees that the integral is a function of just the 
invariant £ = (x — x') 2 /Ayy' . This integral can be simplified by integrating over x. Hence, 

considering the sequence of transformations F{£) — > F{p) — > F{k) defined in (4.18) and 
(4.25), we obtain 



F(sinh 2 9) = fdke- lk9 P 1 (k)F 2 (k). 

27T J 



Ay ' Z (D.2) 



If F{p) can be determined then F{£) can be recovered by using (4.19). For our purposes 
it is necessary to extend this method to deal with integrals which transform as conformal 
tensors. For illustration we first consider 

F®(XpX v -±6 lw )=f™dz Jd^rj^^X^-^F^F.ii'), {D.3) 
_ i 

for X^ = y[£(l + £)] 2 d^£, a conformal vector of scale at x. Again the form of the 
integral is dictated by conformal invariance to be given in terms of the single function 
F(£). To reduce (D.3) to the previous case we introduce the differential operator 



= d^d u + -{n»d v + nvdp) - -8^ (d 2 + -n-d) , {DA) 
which is constructed to give 

If we now set 



F(£) = 4^(1 + (D.Q) 
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then to evaluate (D.3) it is sufficient to calculate 

where T\ is given in terms of F\ similarly to (D.6). From (D.6) and (4.18) we can compute 
the transform T without explicitly finding T by solving (D.6) since 

f-oo 

,A+1 spll i 



T{9) = f(A + 2) J dU U * {U + P) • {D - 8) 

For integrals involving the energy momentum tensor we are interested in integrals for which 
F\ — ► F T where 

F T (0=[^ + Or hd - (D.9) 
In this case the required transforms are particularly simple 

^ (sinh ^) = ^_L_e-<«)ir h(k) = \ Sd wT ^- ? . (D.10) 

It remains to treat integrals of the form given by (5.26) for which we consider the 
following expression 

dz J d^rj-y^X^ - -V) (x' a X' p - -S ap )F T (i)H(i') . (D.ll) 

We now write using (D.5) 

Q^o P = ^yy , fv^V' ap G{i) (D.12) 

where 

Q{i) = So dZ I ^ Wf^^ ' 4 ^ (1 + = H{0 ' (jD ' 13) 
Applying the transformations as described earlier we easily find 

£(sinh 2 9) = — [dk e~ M f T {k)H(k) . {D.IA) 
2-k J 

The result in (D.10) for Pr{k) then allows us to write 

{{d + l) 2 - ^)^(sinh 2 9) = ^S d H(smh 2 9) , (D.15) 
which in turn translates into 

+ 0^ + d« + 1)^ - = -^S d H(Z) ■ (D.16) 
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This equation is more conveniently written as 

(fd + + C + *>« + 1)£ + + 2 )6"(f) = -3*16^+0^ ' (B ' 17) 

To evaluate (D.13) we use (2.17) and 



1 4£ + l ^ v 1 

to obtain 



(D.18) 



£p^<rp = X Mi y,o- p (s) c(£) + (Xp,X v - (^'o^p - ^5ap)a(£) 

+ ^^/^(s) + /x <-> i/, a <-> p (D.19) 
4 , , 4 4 \ 

where the three invariant functions are given by 



c(o = 88" ($ , ho = -8(i + o (i + t±) <?"(o , 



a(0 = 16(1 + 2 (2 + 4£^ + £ 2 ^)s"(£) 



(D.20) 



An important test on the results (D. 19,20) is provided by the conservation equation which 
may be obtained by using (2.25) (with the arbitrary constant c = 0) giving 

9 »^tj)d G ^ ap = ~ d ~d^ Sd (2y) d + 1 j{~( Iu<T ^ X ' p + Iv P^ X ' a ~ \ X v 8 °p) H (£) 

+ 2X V (x' a X' p - - d 8 ap ) {H(0 + £(1 + 0H'(0) 



(£•21) 

This gives two differential equations relating a, 6, c in (D.19) to H which, with the results 
(D.20), are equivalent to (D.17). 

For presenting the results it is more convenient to define from (2.37,38) 
C(0 = - HO ~ 6(0 = 4(1 + 2£)S"(0 + 8(1 + £)£^S"(£) , 



A ® = h [d ~ l) ~ 1)(a(0 + m)) + 



= - l(d - i) 2 (i +oe((2c + 1)^ + 2)^(0 

1 



(D.22) 



46 



For our applications we consider for H(£) functions of the form 

m) = WTW (B ' 23) 

From the definition of Ti{£) in terms of H{£) in (D.13) and using (D.8) we find 

With this result and (D.10) it is sufficient to find the inverse transform of functions of the 
form 

~9ab( k ) ~ r(a- ik)T(a+ \k) 

±k 2 + y? F(b- \k)T{b + \k) 

where we are ultimately interested in taking \i = j(d + 1), a = a — |A and 6 = ^(A + 3). 
By contour integration it is straightforward to obtain 



J(sinh 2 6>) = — / dke- ike ~i{k) 
2tv J 

4r(2a) 



r(6-a)r(6 + a)(^ 2 - a 2 ) 

v V" 1 (2a)n(l + a ~ &)n(a - l*)n{a + e -4(a+n)|fl| i D - 2( >) 

n\ (a + b) n (l + a- fi) n (l + a + fj,) n 



n=0 

r( -//)r(o + ^) 2 M 

T(b-fj,)T(b + n) n 

For the final inverse transform, we observe that it is only necessary to find /"(£) which is 
given by inverting (D.8) 

' ' - -A-3 



This simplifies the calculation significantly. From (B.17) and (B.18) it follows that 
„ (n = _J r(2a + A + 2) 1 v h N (-1) N 

7r A 4 2a-l F ( 6 _ a ) r ( 6 + a )( At 2 _ a 2) ^ + £)]a+*A+l A/ AT! [4£(1 + 

1 r(a-^)r(a + ^) T(2^ + A + 2) 

+ ^42^-1 F ( 6 _ ^ r ( 6 + ^ f^/x + 1) 

X ^7^PTI F ^ + |A + 1, /, - |A - i; 2^ + 1; -^^y) , (D . 28) 

(a - ±A- |)jy(a+ §A + 1) N 
(l + 2a) N 

El (2a) n (l + a - b) n (a - /j) n (a + [i) n (a + §A + §)»(-#)« 



A" 



n= 



£ n! (a + 6) n (l + a - ^) n (l + a + y) n {a |)„(1 + 2a + iV) n ' 
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Now if we use the fact that b = |(A + 3) then the series in simplifies to a terminating 
5 F 4 ( . ; . ; 1) series and may be summed exactly once again by Dougall's theorem [27] giving 

(l + a - /i)jv(l + a + m)at 

With this result the series in (D.28) can be written exactly as a 3-F2 hypergeometric 
function. Putting it all together we obtain for H given by (D.23) with general a 

Q"(0 S<l ' ' 



32d(d-2a)(± + a) [£(1 + £)] 1+Q 
x | 3 F 2 (1, 1 + a, a - ±d; § + a, 1 + a - §d; - — L- ^ ) (D.30) 

r(i+ id-g)r(i + a-id)r(i + |d)r(| + «) r r -id 
r(i + a)r(| + Jd) L ^ 1 + ^ J 

Note that when a = \d the two terms cancel so that there is no pole. Also from (B.28) 
Q"{i) has a leading behaviour for £ — > with terms oc £, £ 1+a . We have checked that (D.30) 
satisfies (D.17) with the last term representing a solution of the homogeneous equation. 

In order to evaluate the integral in (5.26) we need to take, for the ordinary case 
corresponding to (5.29a), 

N 4dA A , ord 1 

m - d -s!Jd^W [ai + 0] ^- (IX31) 

Applying (D.30) gives then 

N d-2 1 



£"(£)ord = 



, d-1 



x {3^2 (1, d - 1, - 2; d - |, \d - 1; (^- 32 ) 

r(id-i)r(3-id)r(id + i)r(d-i) i d _ 2 
r(d-i)r(id+l) ^li + «J 

For C, .4. defined as in (D.22) we have 
N d-2 2£ + l 



c(0 



ord 



I d-1 



^(^-i) 2 [e(i+e)]' 

x {3F2 (1, d - 1, \d - 2; d - |, \d - 1; + (^ 33 ) 



r(id-i)r(3-id)r(id + i)r(d-|) * d _ 2 

r(d-i)r(id+i) + 
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and 



N 8 2 1 

AT 2 1 



,d-l 



K(i + 0]' 

x {3^2 (1, d - 1, \d - 2; d - |, |d; ~ 4|(1 1 +0 ) 

_ r(id) 2 r(3-id)r(d-f) ^ / 

r(d-i)r(id-i) L44U + 4jJ v 



(DM) 



1 + 



d 



d-14C(l + 



Appendix E 

The crucial step in deriving (7.20) is to obtain, with D ap (y,d) defined as in (7.16), 
the result for arbitrary H(£) 

= (d 2 - 1) (Ayy') d \ X ^P^ 2^ + 1) 

+ \ fax' a I„ P (s) +v~v,a~p- - d 5 ap X^X v ) (j^Hit) - H'(0) 

+ X P X V (x' a X' p - - d 8 ap ) fa + 1)#"(£) + 2H'(0 - j^H(O) 

- X P X V (n a n p - - d 5 ap ) fa + 1)H"(£) + d(£ + \)H\t) 

The terms dependent on n a n p violate conformal invariance but they disappear on removing 
the trace in \iv and H(£) — > + l)]~5 d as appropriate for application in (7.19). 
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